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SECTION  1 


INTRODUCTION 

In  the  early  days  of  gas  turbine  development,  the  view  was  widely  held 
that  within  a  few  years  engines  would  be  capable  of  operating  efficiently 
on  a  wide  variety  of  cheap  fuels.  Although  Whittle  had  chosen  kerosine  on 
the  grounds  of  cost,  availability  and  ease  of  handling,  it  was  generally 
believed  that  normal  development  procedures  would  yield  engines  capable  of 
burning  a  wide  range  of  gaseous,  liquid  and  even  pulverized  solid  fuels. 
However,  this  early  optimism  was  short-lived,  at  least  for  aircraft  engines, 
as  experience  soon  revealed  the  important  and  restrictive  effects  of 
aircraft  and  engine  requirements  on  the  physical  and  chemical  properties 
of  the  fuel . 

In  due  course, the  extensive  research  initiated  by  the  major  fuel 
and  engine  companies  and  the  various  governmental  research  establishments 
yielded  increasing  knowledge  of  fuel  properties  and  the  manner  in  which 
they  affect  aircraft  performance,  safety  and  reliability.  This  improved 
understanding  helped  to  eliminate  or  ease  many  formidable  problems,  but 
from  time  to  time  various  new  problems  emerged,  mainly  as  a  result  of  the 
continuing  demand  for  increases  in  aircraft  speed,  range  and  operating 
altitude.  In  the  late  1960s  the  growth  of  public  concern  over  air  pollution 
caused  by  combustion  processes  eventually  culminated  in  legislation  for 
closer  control  of  exhaust  emissions,  notably  in  terms  of  reductions  in 
carbon  monoxide,  unburned  hydrocarbons,  smoke  and  nicric  oxides.  The 
restraining  effects  of  these  requirements,  coupled  with  the  fact  that  the 
aircraft  user  has  traditionally  been  able  to  draw  his  fuel  supplies  from 
the  highest  quality  feedstock,  fostered  conservatism  to  the  extent  that 
current  airline  fuels  do  not  differ  markedly  from  the  kerosine  used  by 
Whittle. 


In  general,  for  any  given  aircraft  application,  the  optimum  fuel  is 
one  that  represents  the  best  compromise  solution  to  the  various  problems 
confronting  the  fuel  companies,  the  engine  and  aircraft  manufacturers  and 
the  operator.  For  civil  aircraft  the  main  requirements  are  safety, 
reliability,  low  cost  and  ease  of  handling.  For  military  aircraft  fuel 
cost  is  of  secondary  importance  compared  with  availability,  supply  logis¬ 
tics  and  the  need  for  trouble-free  operation  over  a  wide  range  of  conditions. 

The  most  dominant  fuel  issues  of  today  are  those  of  cost  and  avail¬ 
ability.  The  steps  now  being  taken  to  ensure  future  supplies  of  fuels  for 
gas  turbines,  in  addition  to  various  measures  of  fuel  conservation,  include 
the  exploitation  of  alternative  fuel  sources  and  the  acceptance  of  a 
broader  specification  for  aviation  fuels.  It  is  clearly  of  paramount 
importance  that  prediction  techniques  be  established  for  estimating 
accurately,  for  any  given  combustor,  the  impact  of  any  change  in  fuel 
specification  on  hardware  durability  and  the  key  aspects  of  combustion 
performance. 

A  complicating  factor  in  the  attainment  of  this  goal  is  that  the 

effect  of  a  change  in  fuel  properties  is  not  constant  for  all  combustors 

but  varies  between  one  combustor  and  another,  due  to  differences  in  operating 

conditions  and  differences  in  design.  For  example,  the  effect  of  an 

increase  in  carbon/hydrogen  ratio  on  liner  wall  temperature  is  much  greater 

for  combustors  featuring  fuel-rich  primary  zones  than  for  combustors  in 

which  the  primary  zone  is  fuel -weak.  This  is  because  with  rich  primary 

zones  most  of  the  heat  transferred  to  the  liner  wall  is  by  radiation, 
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which  is  proportional  to  cT^  .  Thus  liner  wall  temperature  is  dependent 
on  the  flame  emissivity,  e,  which,  in  turn,  is  dependent  on  the  C/H  ratio 
of  the  fuel.  With  fuel -weak  primary  zones,  "however,  most  of  the  heat 


transferred  to  the  liner  wall  is  by  forced  convection.  Here  the  dominant 
term  is  the  gas  temperature,  T  ,  which  is  fairly  insensitive  to  changes  in 
C/H  ratio.  In  consequence,  quite  large  changes  in  C/H  ratio  produce  only 
a  slight  effect  on  liner  wall  temperature. 

Another  complicating  factor  is  that  the  various  properties  and 
characteristics  of  petroleum  fuels  are  so  closely  interrelated  that  it 
is  virtually  impossible  to  change  any  one  property  without  affecting  many 
others.  Thus  the  classical  approach  to  experimental  research,  which  is 
based  on  examination  of  the  effects  of  varying  one  independent  parameter, 
while  maintaining  the  others  constant,  is  precluded  from  the  outset. 

Fortunately,  there  are  several  mitigating  factors  that  help  to  ease 
the  situation.  For  example,  atomization  quality  is  influenced  only  by  the 
physical  properties  of  the  fuel,  namely,  viscosity,  surface  tension,  and 
density,  all  of  which  are  easily  measured  by  standard  laboratory  techniques. 
Moreover,  it  can  also  be  shown  that  evaporation  rates  are  closely  linked  to 
the  physical  properties  of  the  fuel,  for  example,  pp  provides  a  useful 
indication  of  fuel  volatility. 

As  far  as  the  actual  combustion  process  is  concerned,  it  is  found  that 
chemical  reaction  rates  vary  only  slightly  between  the  various  hydrocarbon 
fuels  of  interest  to  the  aircraft  gas  turbine.  This  is  partly  because 
these  fuels  exhibit  only  slight  differences  in  adiabatic  flame  temperature, 
but  also  because  before  entering  the  true  reaction  zone  all  the  fuels  are 
largely  pyrolyzed  to  methane,  other  1-2  carbon  atom  hydrocarbons,  and 
hydrogen.  Hence  the  gas  composition  in  the  reaction  zone  is  substantially 
independent  of  the  parent  fuel.  Thus,  provided  the  discussion  is  restricted 
to  the  anticipated  range  of  aircraft  fuels,  as  reflected  in  Table  1,  any 


differences  that  occur  in  ignition  performance*  lean-blowout  limits  and 
combustion  efficiency,  will  be  caused  mainly  by  differences  in  the  physical 
properties  of  the  fuel  insofar  as  they  control  the  quality  of  atomization 
and  the  ensuing  rate  of  evaporation. 

These  same  physical  properties,  along  with  the  critical  liner  design 
features,  and  the  combustor  operating  conditions,  also  play  a  role  in 
determining  the  rates  of  formation  of  carbon  monoxide  (CO),  unburned 
hydrocarbons  (UHC)  and,  to  a  lesser  extent  N0X-  Smoke,  and  the  soot 
formation  in  the  combustion  zone  that  gives  rise  to  smoke,  are  strongly 
dependent  on  fuel  chemistry.  Thus  fuels  of  different  chemical  composition 
(especially  in  regard  to  hydrogen  and  aromatics  content)  are  found  to 
exhibit  wide  variations,  both  in  radiant  heat  flux  to  the  liner  walls  and 
in  the  level  of  smoke  concentration  in  the  engine  exhaust  gases.  It  follows 
that  the  success  of  any  proposed  relationships  for  luminous  flame  radiation, 
liner  wall  temperature,  and  soot  and  particulate  emissions  will  depend  upon 
accurate  identification  of  the  initial  fuel  components  that  govern  the 
nature  and  rate  of  soot  formation  in  the  flame. 

In  subsequent  sections  the  main  combustor  performance  parameters  are 
discussed  in  turn.  In  each  case  an  outline  is  given  of  the  general  approach 
employed  in  identifying  the  basic  relationships  between  the  relevant  fuel 
properties  and  each  individual  aspect  of  combustor  performance.  Meaningful 
relationships  are  sought,  not  on  the  basis  of  statistical  techniques,  but 
from  an  understanding  of  the  fundamental  physical  and  chemical  processes 
involved.  The  general  approach  has  been  either  to  enhance  existing 
correlations  or  to  replace  them  with  new  correlations  that  are  based  on  a 
firmer  scientific  footing.  It  is  hoped  that  the  relationships  developed 
in  this  program  will  make  a  useful  contribution  to  future  combustor  designs. 
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SECTION  II 


BASIC  DATA 

In  recent  years  the  USAF,  Army,  Navy,  and  NASA,  along  with  engine 
manufacturers , have  initiated  programs  to  determine  the  effects  of 
anticipated  future  fuels  on  existing  engines.  As  a  result  of  these  studies 
data  have  become  available  that  yield  new  and  useful  insights  into  fuel 
property  effects  on  combustion  performance.  These  data,  that  are  contained 
in  references  1  to  6,  provide  the  basic  material  for  this  investigation. 

In  addition  to  a  considerable  body  of  evidence  on  the  effects  of 
fuel  property  variations  on  the  performance,  exhaust  emissions,  and 
durability  characteristics  of  the  combustors  investigated,  references  1 
to  6  also  contain  detailed  information  on  all  the  relevant  chemical  and 
physical  properties  of  the  fuels  employed.  These  fuels  were  supplied  by 
the  USAF  for  combustion  system  evaluation.  They  included  a  current  JP4, 
a  current  JP8,  five  blends  of  the  JP4 ,  five  blends  of  the  JP8  and,  in 
some  cases,  a  No.  2  diesel  fuel.  The  blends  were  intended  to  achieve 
three  different  levels  of  hydrogen  content;  i.e.  12,  13  and  14  percent 
by  mass. 

The  rationale  for  the  diesel  fuel  was  to  approximate  the  Experimental 
Referee  Broad  Specification  (ERBS)  aviation  fuel  that  emanated  from  the 
NASA-lewis  Workshop  on  Jet  Aircraft  Hydrocarbon  Fuel  Technology  [7J. 

The  JP4,  JP8  fuels  and  their  blends  were  chosen  to  span  systematically  the 
possible  fuel  variations  in  key  properties  that  might  be  dictated  in  the 
future  on  grounds  of  availability  and  cost,  and  the  use  of  non-petroleum 
sources  for  jet  fuel  production. 

The  key  chemical  and  physical  properties  of  the  fuels  selected  are 
lisced  in  Tables  1  and  2.  Additional  information  on  the  distillation 
characteristics  of  the  test  fuels  is  contained  in  Fig.  1. 


Properties.  (Reference  6) 
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Figure  1.  Distillation  Characterist 


SECTION  III 


FUEL  ATOMIZATION 


The  quality  of  the  experimental  data  contained  in  references  1  to  6 
is  generally  high.  Although  detailed  information  on  the  main  liner 
dimensions  and  airflow  distribution  is  somewhat  sparse,  it  was  usually 
possible  to  deduce  these  parameters  to  an  acceptable  level  of  accuracy. 
Only  in  one  area,  namely  that  of  fuel  atomization,  did  lack  of  accurate 
information  prove  a  real  impediment  to  the  investigation.  It  is 
strongly  advised  that  in  future  experimental  studies  on  fuel  effects 
every  effort  should  be  made  to  determine  mean  drop  size  and  drop  size 
distribution  for  all  fuels  over  wide  ranges  of  combustor  operating 
conditions.  Sufficient  measurements  should  be  undertaken  to  allow  the 
spray  characteristics  at  intermediate  test  conditions  to  be  obtained  by 
interpolation. 

In  the  absence  of  actual  measured  values,  the  mean  drop  size  (SHD) 
was  calculated  using  one  of  the  following  two  expressions: 

For  airblast  atomizers  [8] 
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For  pressure  swirl  atomizers  [9] 
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Equation  (1)  takes  full  account  of  variations  in  fuel  properties 
(op,  pp  and  Up),  air  properties  (pft  and  \i^)  and  atomizer  geometry  (Dp  and 
Dh).  The  values  of  the  constants  and  exponents  in  this  equation  were 
established  in  a  number  of  experimental  studies  that  covered  much  wider 
ranges  of  fuel  and  air  properties  than  are  needed  for  the  present  investi¬ 
gation.  Thus  the  main  source  of  error  in  the  use  of  Eq.  (1)  stems  from 
uncertainties  surrounding  the  values  to  be  assigned  to  the  atomizer 
dimensions.  With  Eq.  (2)  problems  arise  in  the  calculation  of  pg,  since 
the  primary-zone  temperature  cannot  be  estimated  accurately.  Another 
potentially  serious  source  of  error  that  applies  to  both  equations,  is 
that  all  the  experimental  work  involved  in  their  formulation  was  carried 
out  under  cold,  i.e.  non-burning,  and  fairly  quiescent  conditions. 

Clearly  drop  sizes  could  be  appreciably  different  in  the  true  combustor 
environment  due  to  the  combined  effects  of  high  temperature,  high  turbu¬ 
lence,  and  strong  airflow  currents. 


SECTION  IV 


ANALYSIS 

The  objective  of  the  analysis  was  to  develop  appropriate  relationships 
and  correlations  between  combustion  performance  and  the  relevant  fuel 
properties,  combustor  design  features,  and  combustor  operating  conditions. 
The  method  followed  was  to  study  each  aspect  of  combustion  performance 
from  as  fundamental  a  viewpoint  as  possible,  making  full  use  of  existing 
knowledge  on  the  basic  chemical  and  physical  processes  involved.  Only 
under  circumstances  where  strict  adherence  to  this  dictum  would  clearly 
yield  an  impractical  solution  was  a  less  rigorous  approach  adopted. 

It  was  found  convenient  to  divide  the  analytical  work  into  the 
following  separate  tasks: 

1.  Combustion  efficiency 

2.  Lean  blowout  limits 

3.  Ignition  performance 

4.  Liner  wall  temperature 

5.  Emissions:  carbon  monoxide,  unburned  hydrocarbons,  oxides  of 
nitrogen,  and  smoke 

6.  Pattern  factor 

As  stated  in  the  previous  section,  the  paucity  of  available  data 
on  the  fuel  spray  characteristics  of  the  combustors  described  in  references 
1  to  6  proved  to  be  a  serious  obstacle  to  the  successful  accomplishment  of 
these  tasks.  This  will  become  apparent  in  subsequent  sections  of  the 
report,  where  it  can  be  seen  that  the  quality  of  the  correlations  achieved 
is  distinctly  higher  for  those  aspects  of  performance  where  the  dependence 
on  mean  drop  size  is  either  small  or  ncn-existent.  For  some  performance 
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parameters,  such  as  lean  blowout  and  lean  lightoff  limits,  it  is  possible 
to  express  the  mean  drop  size  for  any  given  fuel  in  terms  of  its  value 
relative  to  that  obtained  with  JP4  for  the  same  fuel  nozzle  at  the  same 
flow  rate,  so  that  any  errors  incurred  through  lack  of  detailed  information 
on  mean  drop  size  are  thereby  greatly  diminished.  However,  for  certain 
other  performance  aspects,  such  as  pattern  factor  and  carbon  monoxide 
emissions,  it  is  important  to  know  for  each  fuel  at  all  the  key  operating 
conditions  the  absolute  value  of  mean  drop  size  in  order  to  calculate 
the  volume  expended  within  the  combustion  zone  in  fuel  evaporation.  In 
the  absence  of  accurate  values  of  SMD  and,  bearir^  in  mind  that  one  of 
the  recommendations  for  the  second  phase  of  this  program  is  that 
measurements  of  SMD  should  be  carried  out  on  all  the  fuel  nozzles  of 
relevance  to  the  investigation  [see  references  1  to  6],  it  was  decided 
that  for  Task  1  the  analysis  of  those  performance  aspects  which  demanded 
a  precise  knowledge  of  SMD  would  be  confined  to  a  relatively  small  number 
of  combustors,  just  sufficient  to  demonstrate  the  validity  of  the 
analytical  procedures  employed.  The  combustors  selected  for  more  detailed 
study  were  the  J79-17A,  F101 ,  TF41 ,  and  J79-17C.  These  combustors  were 
chosen,  partly  because  their  atomization  characteristics  were  fairly  well 
known,  but  also  because  the  relevant  reports  [references  1  to  4]  were 
available  at  the  outset  of  Task  1.  More  detailed  analysis  of  the  TF39, 

J85,  TF33  and  F100  combustors  will  be  performed  as  part  of  Task  II  of  the 
WPAFB  "Fuel  Effects  on  Gas  Turbine  Performance"  program,  which  also 
includes  plans  for  acquiring  a  comprehensive  body  of  drop  size  data  for 
these  combustors. 

The  analyses  that  were  conducted  to  obtain  correlations  for  combustion 
efficiency,  lean  blowout  limits,  lean  lightup  limits,  liner  wall  temperature, 
pollutant  emissions,  and  pattern  factor,  are  described  in  the  following  sections. 
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SECTION  V 


COMBUSTION  EFFICIENCY 


The  main  factors  affecting  the  level  of  combustion  efficiency  are 
evaporation  rates,  mixing  rates,  chemical  reaction  rates  and  the  air 
loading  on  the  combustor.  Thus  combustion  efficiency  may  be  expressed 
as  [10] 
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In  practical  combustion  systems,  the  maximum  rate  of  heat  release  at 
any  given  operating  condition  may  be  governed  either  by  evaporation,  mixing 
or  chemical  reaction,  but  rarely  by  all  three  at  the  same  time.  However, 
over  the  range  of  operating  conditions  where  the  combustion  process  is  in 
transition  from  one  regime  to  another,  it  is  inevitable  that  two  of  the 
three  key  steps  will  participate  in  determining  the  overall  combustion 
efficiency.  Before  exploring  this  situation,  it  is  of  interest  to  examine 
the  separate  effects  on  combustion  efficiency  of  chemical  reaction,  mixing 
and  evaporation. 

1 .  Reaction  Rate-Controlled  Systems 

If  evaporation  and  mixing  rates  are  both  infinitely  fast,  then  Eq.  (3) 
leads  directly  to  the  well-known  e  parameter 
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This  equation  is  ideally  suited,  and  has  been  widely  used,  in  the 
correlation  of  experimental  data  on  combustion  efficiency  obtained  with 
well  atomized,  light  distillate  fuels,  such  as  JP4  and  Jet  A. 

As  the  length  of  the  combustion  zone  is  usually  proportional  to  its 
diameter  or  height,  H,  the  above  equation  may  be  expressed  in  terms  of 
combustor  volume,  Vc,  with  little  loss  of  accuracy,  i.e. 


/ 


P31-75  Vc  exp  (T3/3OO) 
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(5) 


2.  Mixing  Rate-Controlled  Systems 

If  evaporation  and  chemical  reaction  rates  are  both  infinitely  fast, 
Eq.  (3)  becomes 


=  ^(mixing  rate/air  flow  rate) 


Now  the  rate  of  mixing  between  a  turbulent  air  jet  and  the 
surrounding  gas  is  given  by  the  product  of  the  eddy  diffusivity,  the  mixing 
area  and  the  density  gradient.  If  it  is  assumed  that  the  eddy  diffusivity 
is  proportional  to  the  product  of  a  mixing  length,  1,  and  the  turbulent 
velocity  in  the  air  jet,  then 


mixing  rate  (eddy  diffusivity)  x  (mixing  area)  x  (density  gradient) 
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Substituting  in  Eq.  (6)  for 
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yields 
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Mixing  rate 


(7) 


Under  conditions  where  mixing  is  limiting  to  performance,  combustion 
efficiency  will  depend  on  the  ratio  of  the  mixing  rate  to  the  air  flow  rate 
Thus,  by  assuming  that  turbulence  scale  is  proportional  to  combustor  size, 
Eq.  (7)  becomes 


(8) 


3.  Evaporation  Rate-Controlled  Systems 

The  third  case  to  consider  is  when  mixing  and  reaction  rates  are  fast 
enough  for  fuel  evaporation  to  be  the  rate-controlling  step.  Now,  for  a 
volume  of  air,  V,  containing  fuel  drops  whose  initial  drop  size  is  DQ,  the 
average  rate  of  fuel  evaporation  is  given  by  [9] 

mF  =  8(VPF)(k/cp)g(VDo)c>c  ln(1  +  8)0  +  °'22  ReD  °'5)  (9) 


It  is  assumed  that  as  the  fuel  evaporates,  it  instantly  mixes  and 
burns  with  the  surrounding  air.  Thus,  combustion  efficiency  is  obtained 
as  the  ratio  of  the  mass  of  fuel  evaporated  within  the  combustion  zone  to 
the  mass  of  fuel  supplied,  i.e. 
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Equations  (12)  and  (14)  relate  combustion  efficiency  to  combustor 

dimensions  (via  Vc),  combustor  operating  conditions  (via  k  ,  c  and  p  ), 

9  Pg  g 

fuel  nozzle  characteristics  (via  P0)  and  fuel  type  (via  pp  and  B,  or  , 

Values  of  Agfp  are  shown  plotted  in  Figs.  2,  3  and  4.  These  figures 
contain  plots  cf  *eff  versus  T^n  at  three  levels  of  pressure,  namely  100, 

1000  and  2000  kPa,  and  three  levels  of  ambient  temperature,  namely  500, 

1200  and  2000K.  For  each  value  of  temperature  several  lines  are  drawn  to 

represent  different  values  of  UDQ,  where  U  is  the  relative  velocity 

between  the  fuel  drop  and  the  surrounding  gas,  and  DQ  's  the  initial  drop 

diameter.  The  parameter  T^n,  fuel  boiling  temperature  (K)  at  normal 

atmosphere  pressure,  was  chosen  to  represent  the  fuels'  propensity 

for  evaporation.  It  is  recognized  that  no  single  chemical  or  physical 

property  is  completely  satisfactory  for  this  purpose.  However,  the  average  boiling 


point  (50  percent  recovered)  has  much  to  commend  it,  since  it  is  directly 
related  to  fuel  volatility  and  fuel  vapor  pressure.  It  also  has  the 
virtue  of  being  easy  to  measure,  and  is  usually  quoted  in  fuel  specifications. 


From  a  practical  standpoint  the  concept  of  has  considerable 
advantages  since  it  takes  into  account  the  reduced  rate  of  evaporation 
that  occurs  during  the  initial  droplet  heat-up  period,  as  well  as  the 
enhancement  of  fuel  evaporation  rates  due  to  the  effects  of  forced  convec¬ 
tion.  Thus  plots  of  of  the  type  shown  in  Figs.  2  to  4  greatly  simplify 
calculations  on  rates  of  spray  evaporation  and  drop  lifetimes. 

For  any  given  combustor  the  influence  of  fuel  type  on  combustion 
efficiency  can  be  examined  by  defining  a  dimensionless  efficiency  ratio 
which  is  the  ratio  of  the  combustion  efficiency  of  any  alternative  fuel 
'a'  to  that  of  the  baseline  fuel  'b‘,  when  burning  at  the  same  operating 
conditions.  From  Eq.  (12)  we  have,  for  low  ReQ  , 


\  »F  ln(l*Bb) 

<3 


(15) 


This  equation  may  be  used  to  predict  the  change  in  combustion  efficiency 
due  to  a  change  in  fuel  type  provided,  of  course,  that  fuel  evaporation 
is  known  to  be  the  rate-controlling  step.  If  the  level  of  combustion 
efficiency  of  interest  is  high,  say  >  90  percent,  a  more  accurate  prediction 
of  combustion  efficiency  is  obtained  by  rewriting  Eq.  (15)  as 


1-n, 


1  -n. 


%  V 

o F  Da2  lnO+Bb) 
a 


(16) 
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Although  the  derivation  of  Eq.  (16)  from  Eq.  (15)  cannot  be  justified 
on  mathematical  grounds,  the  substitution  of  reciprocal  inefficiency  for 
efficiency  has  become  accepted  practice  because  it  leads  to  a  more 
accurate  result  in  cases  where  combustion  efficiencies  are  close  to 
100  percent. 

For  pressure-swirl  atomizers  the  mean  drop  size  depends  on  the 
surface  tension  and  viscosity  of  the  fuel.  However,  conventional  fuels 
exhibit  merely  slight  differences  in  surface  tension  values  so  that  only 
the  influence  of  viscosity  on  SMD  need  be  considered.  From  Eq.  (2)  we 
have 

SMD  -  MF0*25  (17) 


Substitution  of  SMD  from  Eq.  (17)  into  Eqs.  (15)  and  (16)  respectively 
gives 


nc  PF  uf  C’5  1nO+B,) 
__a  =  rb  a 

\  pf  uF  0,5  ln(l+Bb) 

3  d 


and 


pf  Mr  0-5  ln(l+B  ) 
Pp  uF  °*5  ln(l+Bb) 


(18) 


(19) 


The  practical  utility  of  Eq.  (19)  for  predicting  the  influence  on 
combustion  efficiency  of  a  change  in  fuel  type  is  demonstrated  in  Fig.  5, 
which  contains  experimental  data  obtained  by  Moses  [12]  on  several  fuels, 
using  a  T63  combustor.  By  designating  Jet  A  as  the  baseline  fuel,  the 
combustion  efficiencies  of  the  other  fuels  are  readily  calculated  from 
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Figure  5.  Illustration  of  Use  of  Eq.  (19)  to  Predict  the  Combustion 
Efficiencies  of  Various  Alternative  Fuels  from  Data 
Obtained  with  Baseline  Fuel.  Experimental  Data  from 
Moses  (Ref.  12). 
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£q.  (19}  using  values  for  B,  pf,  and  Up  from  Table  3-  Owing  to  lack  of 
detailed  information  on  the  fuel  spray  characteristics,  the  effect  of 
any  differences  in  mean  drop  size  between  the  various  fuels  cannot  be 
included.  The  results  of  the  calculations  are  shown  in  Fig.  5.  The 
dashed  line  drawn  in  this  figure  describes  the  combustion  efficiency 
of  the  baseline  fuel  as  determined  experimental ly.  The  full  lines 
represent  the  predictions  of  Eq.  (19)  for  the  other  fuels.  The  level 
of  agreement  between  the  predicted  and  measured  values  is  clearly  very 
satisfactory. 

4.  General  Case 

For  the  general  case  in  which  evaporation,  mixing  and  chemical 
reaction  processes  could  all  influence  the  total  rate  of  combustion,  either 
singly  or  in  combination,  an  appropriate  expression  for  combustion 
efficiency  may  be  obtained  by  substituting  Eqs.  (5),  (8)  and  (14)  into 
Eq.  (3)  to  obtain 


P3Vc 

C3fc°o\  .  W30-6  ,  C5 

3l 

Aeff  (aPl/P3)0'5  P30-75  exp  (T3/300) 

The  constant  Cj  in  the  above  equation  is  included  in  order  to  take  into 
account  possible  effects  of  droplet  interaction  on  the  rate  of  fuel 
evaporation.  represents  a  turbulent  diffusion  coefficient,  while 
is  the  collision  factor  for  the  chemical  reaction. 

5.  Reaction  Rate-  and  Evaporation  Rate-Control  led  Systems 

In  practice  it  is  difficult  to  assess  the  importance  of  mixing  rates 
to  combustion  efficiency.  This  is  due  partly  to  lack  of  accurate  information 
on  the  liner  pressure  drop  parameter,  aPl/P3,  but  ^i^y  to  the  fact  that 
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TABLE  3.  Fuel  Properties  used  to  Construct  Fig.  5 


Fuel 

Density 

kg/rrr 

Mass-transfer 

number 

Viscosity 

kg/ms 

Gasoline  (JP  4) 

692 

6.10 

0.00070 

Kerosine  (Jet  A) 

775 

3.75 

0.00129 

Diesel  oil  (OF  2) 

850 

2.80 

0.0030 

I 

Light  fuel  oil 

875 

2.50 

_ 

0.0037 

mixing  rates  tend  to  be  limiting  only  at  conditions  where  the  level  of 
combustion  efficiency  is  close  to  100  percent,  so  that  deficiencies  in 
Derformance  due  to  mixing  are  difficult  to  discern.  For  these  reasons 
it  is  often  preferable  to  express  combustion  efficiency  as  the  product 
of  the  e  efficiency,  nr  »  and  the  evaporation  efficiency,  n  ,  l.e. 

nc  =  n„  x  n  (21) 

c  -o  ce 

The  second  term  on  the  right  hand  side  of  Eq.  (21)  i.e.,  n  from 

ce 

Eq.  (14),  represent:  the  fraction  of  the  fuel  that  is  vaporized  within 
tile  combustion  zone.  For  n  >1,  n.  =  n  ,  and  Eq.  (21)  reverts  to  the 
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0  parameter  which  denotes  the  fraction  of  fuel  vapor  that  is  converted 

into  combustion  products  by  chemical  reaction. 

From  analysis  of  the  available  experimental  data  on  combustion 

efficiency  the  following  expressions  for  n  and  n  were  derived. 

e  ce 


or  =  1  -  exp 
0 


-0.022  P 


1.3 


3 


Vc  exp  (Tc/400) 


fc">A 


(22) 


and 


nr  =  1  -  exp 
ce 


-36x10  P3Vcxeff 


Tc  V  fcmA 


(23) 


It  should  be  noted  in  Eq.  (22)  that  the  pressure  dependence  is  reduced 
from  1.75  to  1.3.  This  is  because  the  exponent  of  1.75  was  determined  from 
tests  conducted  in  the  pressure  range  from  40  to  200  kPa.  In  the  present 
study  the  data  on  combustion  efficiency  were  obtained  at  higher  pressure 
levels,  where  the  intervention  of  mixing  rates  tends  to  lower  the  pressure 
exponent,  in  this  case  to  a  value  of  1.3.  Another  noteworthy  point  is  that 
temperature  dependence  is  now  expressed  in  terms  of  Tc  rather  than  T^. 
Either  may  be  used  (with  suitable  adjustment  to  the  constant)  but  Tc  was 
found  to  give  the  best  correlation.  Tc  is  the  adiabatic  flame  temperature 
in  the  primary  zone,  assuming  complete  combustion  of  the  fuel.  It  is 
calculated  from  the  expression 

Tc  =  T3  +  *Tc 

where  aTc  is  obtained  from  standard  temperature  rise  charts  for  the  fuel 
in  question,  using  appropriate  values  of  P,,  T,  and  q_  (=q  /f  ). 

O  C  QY  £ 
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Values  of  mean  drop  size,  D0  may  be  calculated  from  Eq.  (1)  or  Eq.  {2} 
using  the  information  on  fuel  nozzle  characteristics  supplied  in  references 
1  to  6.  The  very  satisfactory  correlation  of  combustion  efficiency  data 
provided  by  Eq.  (21)  is  demonstrated  in  Figs.  6  to  13.  These  figures 
include  all  the  relevant  data  on  combustion  efficiency  contained  in 
references  1  to  6.  However,  it  should  be  noted  that  for  some  combustors 
the  measurements  of  combustion  efficiency  did  not  include  all  fuels. 
Moreover,  where  both  predicted  and  measured  values  of  combustion  lie 
close  to  100  percent,  the  results  obtained  for  all  fuels  are  virtually 
the  same  and  are  indicated  on  the  graphs  by  a  single  point. 
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Figure  7.  Comparison  of  Measured  and  Predicted  Values  of  Combustion 
efficiency  for  J79-17C  Combustor. 
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Figure  8.  Comparison  of  Measured  and  Predicted  Values  of  Combustion 
Efficiency  for  FI 01  Combustor. 
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Figure  9.  Comparison  of  Measured  and  Predicted  Values  of  Combustion 
Efficiency  for  TF  41  Combustor. 
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Figure  12.  Comparison  of  Measured  and  Predicted  Values  of  Combustion 
Efficiency  for  IF  33  Combustor. 


Combustion  efficiency  (measured). 


98  r  All  Fuels 


A  Idle 
o  Cruise 

(See  Table  B-5,  [6]) 


84  86  88  90  92  94  96  98  100 

Combustion  efficiency  (predicted),  % 


Figure  13.  Comparison  of  Measured  and  Predicted  Values  of  Combustion 
Efnciency  for  FI 00  Combustor. 


SECTION  VI 


LEAN  BLOWOUT 

Until  the  early  1970's  the  problem  of  lean  blowout  was  always  regarded 
as  relatively  minor.  To  a  large  extent  this  was  due  to  the  almost  universal 
use  of  pressure-swirl  spray  atomizers  of  the  duplex  or  dual -orifice  type. 

The  poor  mixing  characteristics  of  these  pressure-swirl  atomizers  create 
many  performance  problems,  not  the  least  being  a  high  rate  of  soot  formation 
in  the  primary  zone.  However,  this  same  poor  mixing  of  fuel  and  primary  air 
has  the  useful  advantage  of  allowing  combustion  to  occur,  albeit  with  low 
combustion  efficiency,  at  mixture  strengths  well  below  the  normal  weak 
extinction  limit.  In  fact,  lean  blowout  limits  in  excess  of  1000  AFR, 
based  on  overall  combustor  values  of  air  and  fuel  flow  rates,  used  to  be 
quite  commonplace.  In  recent  years  the  continuing  trend  toward  improved 
fuel-air  mixing  prior  to  combustion  (for  example,  airblast  atomizers  and 
prechambers)  in  order  to  reduce  the  generation  of  pollutant  emissions, 
especially  N0X  and  smoke,  has  led  to  a  narrowing  of  stability  limits  and  to 
increasing  concern  over  the  attainment  of  satisfactory  lean  blowout  performance. 

For  homogeneous  fuel-air  mixtures  flame  blowout  occurs  when  the  rate 
of  heat  liberation  intheprimary  zone  becomes  insufficient  to  heat  the 
incoming  fresh  mixture  up  to  the  required  reaction  temperature.  With 
heterogeneous  mixtures,  however,  an  additional  factor  is  the  time  required 
for  fuel  evaporation.  For  fuel  sprays  of  low  volatility  and  large  mean 
drop  size  this  time  is  relatively  long  and  is  often  the  main  factor  limiting 
the  overall  rate  of  heat  release.  In  the  analysis  of  lean  blowout  limits 
it  is  convenient  to  consider  homogeneous  mixtures  first  and  then  to  examine 
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how  the  resuits  obtained  should  be  modified  to  take  account  of  fuel 


evaporation. 

For  homogeneous  mixtures  it  has  been  shown  that  the  lean  blowout 
fuel/air  ratio  depends  on  the  inlet  air  velocity,  pressure,  and  temperature, 
and  on  the  size  of  the  combustion  zone  [13].  The  relationship  is  of 
the  form 


(24) 


This  equation  was  used  by  Ballal  and  Lefebvre  [13]  to  correlate 
measured  values  of  weak  extinction  equivalence  ratio  for  propane-air 
flames  stabilized  on  bluff-body  flameholders,  using  experimental  data  to 
determine  optimum  values  of  n,  b  and  x. 

Equation  (24)  may  also  be  used  to  predict  the  lean  blowout  limits  of 
combustion  chambers  suppl ied  with  heterogeneous  fuel-air  mixtures,  provided 
that  the  rate  of  fuel  evaporation  is  sufficiently  high  to  ensure  that  all 
the  fuel  is  fully  vaporized  within  the  primary  combustion  zone.  If  the  fuel  does 
not  fully  vaporize,  then  clearly  the  "effective"  fuel/air  ratio  will  be 
lower  than  the  nominal  value.  However,  if  the  fraction  of  fuel  that  is 
vaporized  is  known,  or  can  be  calculated,  it  can  then  be  combined  with 
Eq.  (24)  to  yield  the  fuel/air  ratio  at  lean  blowout,  i.e.. 


J(heterogeneous) 


'(homogeneous) 


where  is  the  fraction  of  fuel  that  is  vaporized  within  the  primary 
combustion  zone. 
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From  analysis  of  the  factors  governing  the  rate  of  evaporation  of  a 


fuel  spray,  it  was  found  [14]  that 


0.5, 


U  s  8  (p„/pc)  (k/cn)„  ln(l  +  B)  (V  /m  D.  )  (1  +  0.22  Ren  u‘")  (26} 


g  f 


'P  9 


PZ'  pz  0 


8  (k/c  )  in  (1  +  B)  (1  +  0.22  Re,0-5) 
or,  substituting  for  -  - *— * - - - - — - — 2 - 


(27) 


^  _  8  pq  VpzXeff 

Tf  .  o 

™pz°o 


8  pq  Vpz  Aeff 


(28) 


It  should  be  noted  that  Eqs.  (26)  and  (28)  allow  ff  to  exceed  unity. 
When  this  occurs  it  simply  means  that  the  time  required  for  fuel  evaporation 
is  less  than  the  time  available,  so  that  the  fuel  is  fully  vaporized  within 
the  recirculation  zone.  In  these  circumstances  should  be  assigned  a 
value  of  1.0,  so  that 


xBO 


(heterogeneous) 


=  q 


IBO 


(homogeneous) 


(29) 


Substitution  of  q,  an  from  Eq.  (24)  and  f^  from  Eq.  (28)  into 
L80(hom)  f 

Eq.  (25"'  gives 


f  *  O+x)  2 
PZ  mA  Do 

\bo  a  TTPncrrn^ixT 


*pz 


exp(xT3/b)  Aeff 


(30) 


An  additional  term  to  be  incorporated  into  tq.  (30)  is  the  heating 
value  of  the  fuel.  This  stems  from  the  assumption  that  for  any  given 
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operating  condition  the  lean  blowout  temperature  is  the  same  for  all  hydro¬ 
carbon  fuel-air  mixtures.  This  means  that  fuels  with  a  high  gravimetric 
heat  content  should  be  capable  of  burning  down  to  weaker  mixture  strengths 
than  fuels  having  a  lower  heat  content.  With  this  modification  Eq.  (30) 
becomes 


The  first  term  on  the  right-hand  side  of  the  above  expression  is  a 
function  of  combustor  design.  The  second  term  represents  the  combustor 
operating  conditions.  The  third  term  embodies  the  relevant  fuel -dependent 
properties. 

Unfortunately  the  available  experimental  data  do  not  allow  the  values 
of  n  and  x  to  be  determined  with  any  degree  of  accuracy.  This  is  due  mainly 
to  the  fact  that  for  any  given  engine  the  air  mass  flow  rate  is  always 
roughly  proportional  to  air  pressure  which  mokes  it  very  difficult  to 
isolate  their  separate  effects  on  qLBQ.  A"  coat  can  reasonably  be 
deduced  from  Eq.  (31)  is  that  the  exponents  of  m»  and  V  should  be  the 

r\  pZ 

same  and  that  the  pressure  exponent  should  be  somewhat  higher  hy  an  amount 
depending  on  the  effective  reaction  order.  Analysis  of  the  experimental 
data  on  suggests,  in  fact,  that  the  pressure  exponent  is  about  30 
percent  larger  than  that  of  the  air  mass  flow  rate.  Also,  the  temperature 
dependence  corresponds  to  the  relationship 

^LBO  ^  “ ( T^/ 300) 


-  38  - 


Thus  the  simplest  form  in  which  Eq.  (31)  can  be  expressed  that  is 
consistent  with  the  experimental  data  is 


(32) 


where  A"  is  a  constant  whose  value  depends  on  the  geometry  and  mixing 
characteristics  of  the  combustion  zone  and  must  be  arrived  at  experimentally. 
Having  determined  the  value  of  A“at  any  convenient  test  condition,  Eq.  (32) 
may  then  be  used  to  predict  the  lean  blowout  fuel/air  ratio  at  any  other 
operating  condition. 

One  problem  with  Eq.  (32)  is  that  of  assigning  appropriate  values  of 
V  to  all  the  combustors  under  consideration.  Although  values  of  f  , 

r  *•  PZ 

the  fraction  of  the  total  combustor  airflow  entering  the  primary  zone, 
were  either  quoted  directly  in  references  1  to  6,  or  could  easily  be 
estimated,  only  in  one  or  two  cases  could  the  corresponding  primary-zone 
volume  be  assessed  to  the  desired  level  of  accuracy.  To  surmount 
this  problem  it  was  decided  to  substitute  the  predilution  volume,  Vc, 
into  Eq.  (32),  instead  of  This  could  be  justified  on  the  grounds 

that  more  accurate  values  of  Vc  were  available  and,  in  fact,  had  already 
been  used  in  Section  V  for  the  correlation  of  combustion  efficiency  data. 
Moreover,  as  the  ratio  of  primary-zone  volume  to  pre-dilution  volume  tends 
to  be  fairly  constant  for  most  conventional  combustion  chambers,  using 
Vc  instead  of  V has  the  virtue  of  consistency  without  undue  sacrifice 
in  accuracy. 

A  more  serious  drawback  to  Eq.  (32)  is  that  its  solution  demands  an 
accurate  knowledge  of  mean  drop  size  over  the  entire  range  of  operating 
conditions.  As  no  measurements  of  drop  size  were  taken  for  this  investigation, 
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nor  can  they  be  estimated  with  sufficient  accuracy,  this  obstacle  would 
appear  to  rule  out  the  application  of  this  equation  to  the  present  task. 
However,  since  the  main  objective  is  to  study  fuel  effects,  the  problem 
can  be  circumvented  by  rewriting  Eq.  (32)  in  a  different  form,  such  that 
the  fuel -dependent  properties  are  expressed  relative  to  those  of  .;he 
baseline  fuel,  JP4.  With  these  modifications  Eq.  (32)  becomes 

r,  2  L  _  t  )2 


qlB0 


Vc  P3 


A'  fpz  "A 
T73  — 


exp(iy30G) 


LCV„ 


D  at  Tr 


D  at  277.5  K 


(33) 


where  Df  =  mean  drop  size  relative  to  that  for  JP4 

LCV  =  lower  calorific  value  relative  to  that  for  JP4 
r 

\r  =  effective  evaporation  constant  relative  to  that  for  JP4 

at  T  and  P, 
pz  3 

O  *» 

The  term  (D  at  Tp)  / ( D  at  277. 5)c  is  introduced  to  take  into  account  the 
variation  in  drop  size  arising  from  a  change  in  fuel  temperature  from  the 
initial  baseline  value,  which  is  normally  taken  as  277. 5K. 

The  correlations  of  lean  blowout  limits  provided  by  Eq.  (33),  using 
appropriate  values  of  A',  are  illustrated  in  Figs.  14  to  24.  The  close 
agreement  exhibited  between  the  predicted  and  the  measured  values  of  lean 
blowout  fuel/air  ratio  is  clearly  very  satisfactory.  However,  because  the 
measured  values  of  q^g  for  the  F101  seemed  much  too  high,  only  the 
correlation  obtained  with  one  fuel  is  shown  in  Fig.  18  to  demonstrate  the 
ability  of  Eq.  (33)  to  predict  the  effect  of  changes  in  Tg  and  Tp  on  the 
lean  blowout  limits  of  this  combustor.  Also,  it  should  be  noted  that  a 
few  test  points  for  the  F100  are  omitted  from  Fig.  24.  These  points  were 
considered  spurious  since  the  measured  values  of  q^Dg  were  lower  than  the 
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Figure  14.  Comparison  of  Measured  and  Predicted  Values  of  q. Rn 
for  J79-17A  Combustor,  (fuels  i  to  6).  ,jU 
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Figure  16.  Comparison  of  Measured  and  Predicted  Values  of  q 

for  J/9-17C  Combustor,  (Fuels  2A  to  7A) .  LBO 
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Figure  21. 


Comparison  of  Measured  and  Predicted  Values  of  q.ng 
for  TF  39  Combustor.  Plotted  Points  Represent 
Average  of  Five  Values  for  Fuel  and  Air  Temperatures 
between  255  and  277K. 
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24.  Comparison  of  Measured  and  Predicted  Values  of 
for  F 1 00  Combustor. 
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corresponding  predicted  values  (and  other  measured  values  at  similar  test 
conditions)  by  an  order  of  magnitude. 

It  should  also  be  mentioned  that  the  original  correlation  obtained  for 
the  .I79-17A  combustor  was  extremely  poor.  After  further  examination  of  the 
data  it  was  concluded  that  the  designation  of  fuels  in  Tables  C-8,  C-9, 

C-10  and  C-ll  of  reference  1  was  incorrect  and  that  some  mix-up  of  fuels 
had  occurred  during  this  portion  of  the  test  program.  From  inspection  of 
Eq.  (33)  it  is  apparent  that  q^  is  strongly  influenced  by  the  mean  drop 
size  of  the  spray  which,  in  turn,  is  very  dependent  on  fuel  viscosity,  Vp. 
Thus,  to  a  first  approximation,  lean  blowout  limits  should  correlate  with 
fuel  viscosity,  and  this  is  confirmed  by  the  upper  half  of  Fig.  25  which 
demonstrates  an  excellent  correlation  between  these  two  parameters  for 
the  J79-17C  combustor.  The  lower  part  of  this  figure  shows  similar 
plots  for  the  J79-17A  combustor  but,  in  this  instance,  there  is  clearly 
no  correlation  between  and  q Q ,  which  tends  to  suggests  that  the 
reported  fuel  designations  are  incorrect.  To  remedy  this  situation 
the  fuels  were  renumbered,  as  indicated  in  Table  4,  in  a  manner  designed 
to  provide  a  better  correlation  between  Vp  and  ^lbO’  Proce<^ure 

to  a  very  satisfactory  correlation  between  the  experimental  data  on  qLgQ 
and  the  corresponding  predictions  of  Eq.  (33),  as  illustrated  in  Figs. 

14  and  15. 

It  should  be  noted  that  the  observations  made  above  apply  only  to  the 
lean  blowout  data,  and  there  is  no  reason  to  suspect  the  published  fuel 
numbers  for  lean  lightoff  limits  in  Tables  C-8  to  C-ll  of  reference  1. 

The  full  lines  drawn  in  Fig.  26  represent  the  results  of  calculations, 
using  Eq.  (33),  on  the  effect  of  variation  in  initial  fuel  temperature  on 
\bo-  the  experimental  data  show  some  scatter,  inspection  of 
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gure  25.  Relationship  Betwei  Fuel  Viscosity  and  lean  Blowout  limits. 
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Viscosity  at  300  K,  mm2/s 


j  85  Coid  Ground  Start  (20%  rpm) 


the  figure  indicates  reasonably  good  agreement  between  theory  and 
measurement. 

For  each  combustor  a  value  of  A' was  chosen  for  insertion  into 
Eq.  (33)  that  would  provide  the  best  fit  to  the  experimental  data.  These 
values  of  A'  are  given  in  Table  5.  It  was  anticipated  that  similar  types 
of  primary  zones  would  yield  similar  values  of  A',  which  would  clearly  be 
advantageous  from  the  standpoint  of  predicting  the  lean  blowout  limits 
of  future  combustor  designs.  Although  the  variation  in  the  values  of  A' 
exhibited  in  Table  5  virtually  prohibits  such  extrapolation,  it  should 
be  borne  in  mind  that  these  values  embody  all  the  errors  incurred  in  the 
estimates  of  combustion  volume,  the  fraction  of  air  involved  in  primary 
combustion,  and  the  mean  drop  size.  If  these  parameters  could  be 
established  accurately,  then  it  should  be  possible  to  determine  an 
appropriate  value  of  A1  for  each  characteristic  primary  zone/nozzle 
configuration.  By  combining  A' with  f  the  deviation  is  reduced  somewhat, 
as  illustrated  in  Table  5.  The  exceptionally  low  value  of  A‘fpz  for  the 
TF41  combustor  can  be  attributed  to  its  excellent  atomization  character¬ 
istics  at  low  fuel  flows  stemming  from  the  use  of  an  exceptionally  low 
primary  nozzle  flow  number.  Initially,  an  unusually  high  value  of  A'f 
was  calculated  for  the  FI 01  combustor.  This  could  be  the  result  of  an 
error  in  the  recording  of  either  fuel  flow  rate  or  air  flow  rate,  due 
probably  to  the  fact  that  the  tests  were  conducted  on  a  54°  segment  of 
the  full  annular  combustor.  If  the  tabulated  results  in  reference  [2] 
were  accepted  at  face  value  they  would  imply  that  the  engine  was  in  danger 
of  flame-out  at  take  off,  which  is  highly  unlikely.  Dividing  the  published 
values  of  qLBQ  by  (360/54)  reduces  A’f  to  32  which  is  fully  consistent 
with  the  results  obtained  for  the  other  combustors. 
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TABLE  4.  Re-assigned  Fuel  Numbers  for  Lean  Blowout  Data  in 
Tables  t-8  to  C-ll  [1] 


Original  Fuel  Number 

1 

2 

3 

B 

5 

6 

I 

8 

9 

~ — - 

10 

11 

12 

13 

New  Fuel  Number 

1 

12 

10 

9 

8 

11 

2 

5 

3 

6 

13 

1 

1 

TABLE  5.  Values  of  A'  and  B‘ employed  in  Eqs.  (33)  and  (37) 


Engine 

A1 

ftV 

B' 

BV 

J79-17A 

ur> 

CO 

42 

870 

199 

J79-17C 

100 

31 

- 

- 

FI  0 1 

11 

79 

32 

222 

90 

TF  41 

46 

13 

228 

634 

TF  39 

123 

37 

704 

211 

J85 

213 

64 

597 

..  ..  _ 

179 

TF  33 

85 

25.4 

726 

274 

F100 


65 


483 


170 


SECTION  VII 


IGNITION 

1 .  The  Ignition  Process 

For  the  proper  interpretation  of  experimental  data  on  lightup 
performance  it  is  important  to  recognize  that  the  ignition  process  is  not 
a  simple,  single-step  mechanism,  but  actually  occurs  in  two  or  three 
distinct  phases  [9].  Phase  1  is  the  formation  of  a  kernel  of  flame  of 
sufficient  size  and  temperature  to  be  capable  of  propagation.  Phase  2 
is  the  subsequent  propagation  of  flame  from  this  kernel  to  all  parts  of 
the  primary  zone.  Phase  3,  which  applies  only  to  tubular  and  can-annular 
chambers,  is  the  spread  of  flame  from  a  lighted  liner  to  an  adjacent 
unlighted  liner.  Failure  of  any  single  phase  is,  of  course,  equivalent 
to  failure  to  ignite. 

The  success  or  failure  of  Phase  1  depends  mainly  on  whether  or  not 
the  rate  wf  heat  release  by  combustion  in  the  spark  kernel  exceeds  the 
rate  of  heat  loss  to  the  surroundings  by  radiation,  thermal  conduction, 
and  turbulent  diffusion.  The  rate  of  heat  release  is  governed  mainly  by 
the  effective  fuel/air  ratio  adjacent  to  the  igniter  plug,  which  should  be 
close  to  stoichiometric,  and  by  the  size  and  temperature  of  the  kernel, 
which  are  in  turn  determined  by  the  energy  and  duration  of  the  spark.  The 
rate  of  heat  loss  from  the  kernel  is  largely  dictated  by  the  local  conditions 
of  velocity  and  turbulence  and  by  the  quantity  of  excess  fuel  present  in 
the  ignition  zone. 

The  success  of  Phase  2  is  governed  partly  by  the  location  of 
the  igniter,  since  this  determines  whether  or  not  the  hot  kernel  is 
entrained  into  the  primary-zone  reversal  or  is  swept  away  downstream.  It 
is  also  affected  by  all  the  factors  that  control  flame  stability.  Thus, 
an  increase  in  pressure  and/or  temperature,  a  reduction  in  primary  zone 
velocity  ,  or  any  change  in  fuel/air  ratio  toward  the  stoichiometric  value, 

all  of  which  are  beneficial  to  stability,  will  also  improve  phase  2. 
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The  location  of  the  interconnector  is  of  prime  importance  in  phase 
3.  Ideally,  each  tube  entrance  should  coincide  with  the  region  of  highest 
gas  temperature  in  the  liner,  whereas  the  tube  exit  should  be  sited  so  as 
to  ensure  that  the  issuing  hot  gas  flows  directly  into  the  recirculation 
zone  of  the  adjacent  liner.  Phase  3  is  enhanced  by  the  use  of  interconnectors 
in  which  the  flow  area  is  made  large  to  facilitate  the  passage  of  flame, 
and  whose  length  is  kept  short  to  minimize  heat  loss. 

Although  ignition  of  a  combustible  mixture  may  be  accomplished  in  a 
variety  of  ways,  in  the  gas  turbine  it  is  usually  effected  by  means  of  an 
electric  spark,  and  large  amounts  of  energy  are  needed  to  ignite  the 
heterogeneous  and  highly  turbulent  mixture  flowing  at  velocities  of  the 
order  of  25  m/s. 

In  recent  years  a  series  of  detailed  experimental  studies  has  been 
carried  out  on  the  influence  of  electrical  and  flow  parameters  on  minimum 
spark  energy  in  flowing  mixtures  of  fuel  drops  and  air.  These  studies  have 
led  to  a  better  conceptual  understanding  of  the  basic  ignition  process 
and  have  provided  a  useful  theoretical  foundation  for  relating  ignition 
characteristics  to  all  the  operating  variables  involved.  They  confirm 
practical  experience  in  showing  that  ignition  is  made  easier  by  increases 
in  pressure,  temperature  and  spark  energy,  and  is  impeded  by  increases 
in  velocity  and  turbulence  intensity.  They  also  emphasize  that  ignition 
performance  is  markedly  affected  by  fuel  properties  through  the  way  in 
which  they  influence  the  concentration  of  fuel  vapor  in  the  immediate 
vicinity  of  the  igniter  plug.  These  influences  arise  mainly  through  the 
effect  of  volatility  on  evaporation  rates,  but  also  through  the  effects  of 
surface  tension  and  viscosity  on  mean  fuel  drop  size.  The  amount  of  energy 
required  for  ignition  is  very  much  larger  than  the  values  normally 
associated  with  gaseous  fuels  at  stoichiometric  fuel/air  ratio.  Much  of 
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this  extra  energy  is  absorbed  in  the  evaporation  of  fuel  droplets,  the 
actual  amount  depending  on  the  distribution  of  fuel  throughout  the  primary 
zone  and  on  the  quality  of  atomization. 

2.  Theory 

One  model  for  the  ignition  of  fuel  sprays  is  based  on  the  assumption 
that  chemical  reaction  rates  are  infinitely  fast  and  that  the  onset  of 
ignition  is  limited  solely  by  the  rate  of  fuel  evaporation  [15].  Support 
for  this  notion  may  be  found  in  the  literature  on  the  ignition  of  turbojet 
combustors.  For  example,  the  reports  by  Foster  and  Straight  [16]  and 
Wigg  [17]  contain  ample  evidence  that  various  fuel  spray  characteristics, 
such  as  mean  drop  size  and  volatility,  can  appreciably  affect  the  energy 
required  for  ignition.  These  effects  are  due  to  their  influence  on  fuel 
evaporation  rates  that  govern  the  mixture  strength  in  the  ignition  zone. 
Further  confirmation  of  the  importance  of  fuel  vapor  concentration  to 
ignition  is  provided  in  the  more  basic  studies  conducted  by  Rao  and 
Lefebvre  [18]  on  the  ignition  of  heterogeneous,  flowing  kerosine/air 
mixtures.  They  found  that  for  mixtures  weaker  than  stoichiometric  the 
main  Factor  limiting  ignition  is  a  deficiency  of  fuel  vapor  in  the 
ignition  zone. 

The  process  of  ignition  is  envisaged  to  occur  in  the  following 
manner.  Passage  of  the  park  creates  a  small,  roughly  spherical, 
volume  of  air  (the  spav  .<  ernel)  whose  temperature  is  sufficiently  high 
to  initiate  rapid  e.-af  'al,ion  of  the  fuel  drops  contained  within  the 
volume.  It  is  assumed  that  reaction  rates  and  mixing  times  are  infinitely 
fast  so  that  any  fuel  vapor  created  within  the  spark  kernel  is  instantly 
transformed  into  combustion  products  at  the  stoichiometric  flame  tempera¬ 
ture.  If  the  rate  of  heat  release  by  combustion  exceeds  the  rate  of  heat 
loss  at  the  surface  of  the  inflamed  volume,  then  the  spark  kernel  will 
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grow  in  size  to  fill  the  entire  combustion  volume.  If,  however,  the  rate 
of  heat  release  is  less  than  the  rate  of  heat  loss,  the  temperature  within 
the  spark  kernel  will  fall  steadily  until  fuel  evaporation  ceases  altogether. 
This  concept  leads  to  the  definition  of  'quenching  distance'  as  the  critical 
size  that  the  inflamed  volume  must  attain  in  order  to  propagate  unaided, 
while  the  amount  of  energy  required  from  an  external  source  to  attain 
this  critical  size  is  termed  the  "minimum  ignition  energy." 

Analysis  of  the  relevant  heat  transfer  and  evaporation  processes  [19] 
yields  the  following  expression  for  the  quenching  distance  of  quiescent  or 
slow-moving  monodisperse  mists: 


V 


»F»0 


1  n  ( 1  +  3^} 


0.5 


(34) 


It  should  be  noted  that  the  above  equation  was  derived  directly  from 
basic  considerations  of  the  mechanisms  of  heat  generation  within  the  kernel 
and  heat  loss  from  its  surface  and  contains  no  experimental  or  arbitrary 
constants. 

Equation  (34)  and  similar  expressions  for  polydisperse  sprays  of  the 
type  provided  by  most  practical  atomizing  devices,  provide  simple  relation¬ 
ships  between  quenching  distance  and  the  key  spray  properties.  Essentially, 
they  state  that  quenching  distance  is  directly  proportional  to  drop  size 
and  is  inversely  proportional  to  the  square  root  of  gas  pressure. 

An  increase  in  if  and  reduction  in  Pp  both  reduce  dq  because  they  promote 
evaporation  by  increasing  the  surface  area  of  the  fuel.  Similarly,  an 
increase  in  8  will  also  accelerate  evaporation  and  thereby  decrease  dq. 
Values  of  Em^n  may  be  obtained  by  inserting  the  calculated  values  of  d^ 
from  Eq.  (34)  into  the  following  expression 


Emin  =  cpA  °A  ATst  {*/6)  V 


(35) 
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Subsequently,  Ballal  and  tefebvre  [19]  extended  the  model 
described  above  to  include  (1)  the  effects  of  finite  chemical  reaction 
rates,  which  are  known  to  be  significant  for  very  well -atomized  fuels  at 
low  pressures  and  low  equivalence  ratios,  and  (2)  the  effect  of  heat  loss 
from  the  spark  kernel  by  turbulent  diffusion.  Thus  the  model  has 
general  application  to  both  quiescent  and  flowing  mixtures  of  air  with 
either  gaseous,  liquid  or  evaporated  fuel,  or  any  combination  of  these 
fuels. 

From  a  gas  turbine  viewpoint,  the  value  of  the  type  of  model  described 
above  lies  not  so  much  in  its  ability  to  predict  minimum  ignition-energy 
requirements,  since  the  available  spark  energy  is  determined  by  the 
specifications  of  the  igniter  plug  and  the  high-energy  unit;  rather,  its 
importance  is  in  highlighting  the  key  parameters  that  control  ignition  and 
in  providing  quantitative  relationships  or  "rates  of  exchange"  between 
these  parameters. 

Consider,  for  example,  a  combustor  that  normally  operates  on  kerosine 
fuel.  It  is  required  to  know  what  improvement  in  atomization  quality 
would  be  needed  to  achieve  the  same  ignition  performance  when  diesel  oil 
is  burned.  Now  the  relevant  fuel  properties  are  density  and  volatility, 
the  latter  being  represented  in  the  quenching-distance  equations  by  the 
mass-transfer  number  B.  From  Eq,  (34)  we  have 


Thus,  for  unchanged  ignition  performance,  i.e.,  unchanged  d^,  we  can  write 
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(D0)diesel 

(Dc)ker 


[In  (1  +  Bst)/pF]  diesel 


[in  (l  +  Bst)/PF'}  ker 


0.5 


(36) 


Now  for  diesel  oil,  pp  =  900  and  Bst  =  2.80;  for  kerosine,  pp  =  775  and 
B  t  =  3.75.  Substituting  these  values  into  Eq.  (36)  gives 


^o^diesel  _ 
^Vker 


0.86 


which  suggests  that  for  diesel  oil  the  SMD  should  be  reduced  by  about 
20  percent  to  retain  the  ignition  performance  obtained  with  kerosine. 

3.  Data  Analysis 

Application  of  the  theoretical  concepts  described  above  to  the  analysis 
of  the  experimental  data  on  ignition  contained  in  references  1  to  6  was 
inhibited,  as  for  the  "lean  blowout  data,  by  lack  of  accurate  information 
on  mean  drop  size.  So  again  recourse  was  made  to  relationships  in  which 
the  key  fuel  properties  are  not  expressed  in  absolute  terms  but  in  values 
relative  to  those  of  the  baseline  fuel,  JP4.  This  approach  leads  to  the 
following  equation  for  lean  lightup  fuel/air  ratio. 


B'f 


pz 


m. 


D. 


TL0 


P3V-5  exp  (T3/300) 


r 

Tar 


D  at  T r 


D  at  27775K" 


(37) 


This  equation  is  virtually  identical  to  Eq.  (13)  except  for  a 

15  13 

higher  pressure  dependence;  P3  '  versus  P3  . 

The  correlation  of  lightup  data  obtained  with  Eq.  (37)  is  illustrated 
in  Figs.  27  to  38.  The  leve1  of  agreement  between  predicted  and  experimen¬ 
tal  values  is  considered  satisfactory,  especially  in  view  of  the  well  known 
lack  of  consistency  that  usually  characterizes  experimental  data  on  spark 
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Ignition.  The  values  of  B'  and  B'f  for  the  various  combustors  are  listed  in 
Table  5. 

A  curious  anomaly  exists  in  regard  to  the  ignition  data  obtained  from 
the  079  17C  combustor.  For  this  combustor  the  results,  as  shown  plotted 
in  Figs.  29  and  30,  indicate  a  correlating  parameter  of  the  form 


qLL0  =  C7  + 


Vc  P3 


C8  fpz  mA 

1.5 


exp(T3/300) 


V  LCVr 


D  at  Tr 


0  at  277. 5K 


(38) 


Assigning  values  to  C -j  and  Cg  of  3.92  and  200  respectively  provides 
the  excellent  correlation  illustrated  in  Figs.  31  and  32. 

For  the  FI 01  combustor  it  was  found  that,  in  common  with  the  lean 
blowout  data  discussed  earlier,  the  reported  lightoff  fuel/air  ratios 
were. unusually  high.  Thus  only  the  results  obtained  with  the  fuel  having 
the  highest  number  of  test  points  (No.  12)  are  shown  in  Fig.  33  to 
illustrate  the  level  of  prediction  attained.  As  for  the  lean  blowout  case, 
it  was  again  assumed  that  a  computational  error  had  arisen  in  the  calculation 
of  fuel/air  ratios  due  to  tests  being  performed  on  a  54°  sector  rig,  and  a 
correction  factor  of  (54/360)  was  again  used  to  calculate  the  values  of 
B  and  Bf  shown  in  Table  5  for  the  FT 01  combustor. 


(measured),  g/kg 


Fiyut'o  29.  Correlation  of  Experimental  Data  on  q,,n  for  J79-17C 
Combustor .  (fuels  1A  to  6A) . 
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(measured),  g/kg 


(measured),  g/kg 


Figure  3a  Comparison  of  Measured  and  Predicted  Values  of  q 
for  F100  Ccmbestor. 
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SECTION  VIII 
LINER  WALL  TEMPERATURE 


For  the  purpose  of  analysis  a  liner  may  be  regarded  as  a  container 
of  hot  flowing  gases  surrounded  by  a  casing  in  which  air  is  flowing  between 
the  container  and  the  casing.  Broadly,  the  liner  is  heated  by  radiation 
and  convection  from  the  hot  gases  inside  it,  and  is  cooled  by  radiation 
to  the  outer  casing  and  by  convection  to  the  annulus  air.  The  relative 
proportions  of  the  radiation  and  convection  components  depend  upon  the 
geometry  and  operating  conditions  of  the  system.  Under  equilibrium 
conditions  the  liner  temperature  is  such  that  the  internal  and  external 
heat  fluxes  at  any  point  are  just  equal.  Loss  of  heat  by  conduction 
along  the  liner  wall  is  comparatively  small  and  usually  may  be  neglected. 
Under  steady-state  conditions  the  rate  of  heat  transfer  into  the  wall 
must  be  balanced  by  the  rate  of  heat  transfer  out.  Thus,  we  have 

R-j  +  C-|  -  R2  +  1*2  (39) 


1.  Internal  Radiation,  R^ 

This  is  the  component  of  heat  transfer  that  is  most  affected  by  a 
change  in  fuel  type.  It  is  given  by  [20] 


R,  =  0.5  0  (1 


+  C  ) 


T  ^ ^  {7  2.5 

9  g 


-  T  2'5) 
w  ' 


where 


a  =  Stefan  Boltzmann  constant 

e  =  liner  wall  emissivity 
w 

Eg  =  gas  emissivity 
Tq  ~  gas  temperature 
Tw  =  wall  temperature 


(40) 
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The  'bulk'  or  mean  gas  temperature,  Tg,  is  obtained  as  the  sum  of 
the  chamber  entry  temperature,  Tg,  and  the  temperature  rise  due  to  combus- 

ti0"-  SW 

Thus: 

Tg  "  T3  +  &Tcomb 

Alcomb  derived  from  standard  temperature  rise  curves.  The 

appropriate  value  of  fuel/air  ratio  is  the  product  of  the  local  fuel/air 
ratio  and  the  local  level  of  combustion  efficiency.  Most  heat  transfer 
calculations  are  carried  out  at  high  pressure  conditions  where  it  is 
reasonable  to  assume  a  combustion  efficiency  of  100  percent. 

For  the  luminous  flames  associated  with  the  combustion  of  hetero¬ 
geneous  fuel -air  mixtures,  the  value  of  eg  for  insertion  in  Eq.  (40) 
is  obtained  as  [20] 

tg  =  1  -  expC-290  P3  l  (q  t/*5  T^1*5]  (4l) 

where  q  is  the  local  fuel/air  ratio  and  is  the  'beam  length'  of  the 
radiating  gas.  Beam  length  depends  on  the  size  and  shape  of  the  gas 
volume.  For  most  practical  purposes  it  is  given  to  sufficient  accuracy 
[20]  by  the  expression 

s.^  =  ( volume/surface  area) 

The  luminosity  factor,  l,  is  an  empirical  correction  introduced  to 
obtain  reasonable  agreement  between  experimental  data  cn  gas  radiation  and 
predictions  from  Eq.  (40).  Experiments  have  shown  that  luminosity 
factor  depends  largely  on  the  carbon -to -hydrogen  mass  ratio  of  the  fuel 
[9,  20,  21]. 
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The  original  equation  for  L  is  [20] 

L  »  7.53(C/H  -  5.5}0'84  (42) 

Later  the  following  expression  was  suggested  by  Kretschmer  and 
Odgers  [21] 

L  =  0.0691  [C/H  -  1.82]2'71  (43) 

Another  correlation,  which  is  simpler  and  probably  no  less 
accurate,  is  [9] 

L  =  3  (C/H  -  5.2)0'75  (44) 

More  recent  investigations  have  tended  to  emphasize  fuel  hydrogen 
content,  rather  than  carbon/hydrogen  ratio,  as  the  property  most  relevant 
to  flame  radiation.  Figure  39  shows  the  correlation  obtained  by  Blazowski 
and  Jackson  [22,  23]  between  hydrogen  content  and  liner  wall  temperature 
for  several  engines.  The  data  shown  represent  cruise  conditions  with 
combustor  inlet  temperatures  ranging  from  547  to  756  K.  The  parameter 
used  to  correlate  the  experimental  data  is 

(T.  *  T,  )/(T.  -  T3) 

in  which  the  numerator  represents  the  increase  in  maximum  liner  temperature 
over  that  obtained  with  a  baseline  fuel  containing  14.5  percent  hydrogen. 

The  excellent  correlation  of  data  exhibited  in  Fig.  39  could  not  be 

duplicated  for  the  results  contained  in  references  1  to  4,  as  shown  in 

Figs.  40  to  43.  This  is  believed  to  be  due  to  the  fact  that  the  magnitude 

of  the  parameter  (T.  -  T.  )/(T.  -  T3)  is  very  sensitive  to  the  value 

of  T.  .  Since  maximum  values  of  wall  temperature  are  notoriously  difficult 
Lo 
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Figure  40.  Effect  of  Fuel  Hydrogen  Content  on  Liner  Temperature 
Parameter  for  J79-17A  Combustor. 
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Figure  41.  Effect  of  Fuel  Hydrogen  Content  on  Liner  Temperature 
Parameter  for  J79-17C  Combustor. 
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Figure  43.  Effect  of  Fuel  Hydrogen ‘Content  on  Liner  Temperature 
Parameter  for  IF  41  Combustor. 


to  determine,  this  would  appear  to  be  a  serious  drawback  to  the  use  of 
this  type  of  temperature  parameter. 

An  alternative  approach  attempted  here  was  to  find  a  relationship 
between  fuel  hydrogen  content  and  luminosity  factor.  Analysis  of  the 
experimental  data  led  to  the  following  expression 

L  =  336/(percent  hydrogen)*"  (45) 

Substitution  of  this  value  of  L  into  Eq.  (41)  allows  calculations 
of  flame  radiation  to  be  carried  out  for  all  fuels  over  the  entire  range 
of  test  conditions. 

2.  External  Radiation,  R2 

The  radiation  heat  transfer  from  the  liner  wall  to  the  outer  casing, 

R2»  can  be  approximated  by  assuming  gray  surfaces  with  emissivities  ew 
and  e  ,  and  assuming  that  T  and  T  are  approximately  uniform  in  the 
axial  direction.  The  significance  of  R2  increases  with  liner  wall  temperature, 
and  at  low  values  it  can  often  be  neglected.  It  can  be  estimated  only 
approximately  because  of  lack  of  accurate  knowledge  of  wall  emissivities. 

For  this  reason  it  is  sufficient  to  use  the  cooling-air  temperature,  T-^, 
in  place  of  the  unknown  temperature  of  the  outer  air  casing.  Also,  for 
radiation  across  a  long  annular  space,  the  geometric  shape  factor  can  be 
assumed  equal  to  unity.  The  expression  for  the  net  radiation  flux,  then 
reduces  to 


R2  =  o 


ew  £c 


r  +  r  (1 

c  w  ' 


-  c 


75 

c'  A 


<Tw  '  T34> 


(4  6) 
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For  a  tubular  chamber,  (A^/A^  is  equal  to  the  ratio  of  liner  to  casing 
diameter  at  the  section  considered.  For  tubo-annular  systems,  where  the 
depth  of  the  annulus  varies  from  point  to  point  around  the  liner,  an 
average  value  of  0.8  is  used.  For  an  annular  chamber  the  ratio  (Aw/A(.) 
is  slightly  greater  than  unity  for  the  inner  liner  and  slightly  less  than 
unity  for  the  outer  liner. 

Accurate  values  of  emissivity  for  various  materials  may  be  obtained  from 
McAdams  [24].  However,  for  most  practical  purposes  the  following  expression, 
based  on  typical  values  of  emissivity  and  diameter  ratio,  will  suffice: 

R2  =  0.4  0(Tw4  -  T34)  m 

3.  Internal  Convection,  C-j 

Of  the  four  heat  transfer  processes  which  together  determine  the 
liner  temperature,  this  component  is  the  most  difficult  to  estimate 
accurately.  In  the  primary  zone  the  gases  involved  are  at  high  temperature 
and  undergoing  rapid  physical  and  chemical  change.  Further  difficulty 
is  introduced  by  the  existence  within  the  primary  zone  of  steep  gradients 
of  temperature,  velocity  and  composition.  Uncertainties  regarding  the 
airflow  pattern,  the  state  of  the  boundary-layer  development  and  the 
effective  gas  temperature  make  the  choice  of  a  realistic  model  almost 
arbitrary. 

In  the  absence  of  more  exact  data  it  is  reasonable  to  assume  that 
some  form  of  the  classical  heat-transfer  relation  for  straight  pipes  will 
hold  for  conditions  inside  a  liner,  using  a  Reynolds  number  index  consistent 
with  established  practice  for  conditions  of  extreme  turbulence.  This  leads 
to  an  expression  of  the  form  [20] 
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4.  External  Convection,  C2 
This  is  obtained  as  [20] 


C2  =  0.020 
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The  fluid  properties  are  evaluated  at  the  annulus  air  temperature,  Tg. 
In  practice  the  cooling-air  temperature  increases  during  its  passage 
downstream,  but  normally  this  amounts  to  no  more  than  a  few  degrees  and 
can  reasonably  be  neglected. 

For  equil ibrium 


^1  +  ^1  =  +  ^2 


(39) 


Solution  of  this  equation  yields  the  wall  temperature,  T  . 

W 

The  value  of  Tw  as  determined  by  the  method  outlined  above  represents 
the  liner  wall  temperature  that  would  be  obtained  in  the  absence  of  internal 
wall  cooling.  Unfortunately  references  1  to  6  do  not  contain  the  detailed 
information  needed  to  estimate  film-cooling  effects  on  T  .  Thus,  and 
bearing  in  mind  the  lengthy  and  tedious  nature  of  the  procedures  involved, 
it  was  decided  to  calculate  'uncooled'  wall  temperatures  for  four  combustors 
only,  namely  J79-17A,  079-17C,  FlOl  and  T41 ,  in  order  to  ascertain  if 
the  results  obtained  reflected  anticipated  trends  in  regard  to  the  effect 
of  fuel  hydrogen  content  on  liner  wall  temperature.  The  results  of  these 
calculations  are  shown  in  Figs.  44  to  47  for  all  fuels  as  plots  of  Tw  versus 
hydrogen  content. 
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It  may  be  noted  in  Figs.  44  to  4?  that  the  calculated  values  of  7^ 
are  generally  higher  than  the  corresponding  measured  values  due  to  neglect 
of  internal  wall  coding.  Only  at  low  power  conditions,  where  the  errors 
incurred  through  neglect  of  internal  wall  cooling  are  partially  balanced 
by  the  assumption  of  100  percent  combustion  efficiency  in  the  combustion 
zone,  do  the  measured  and  calculated  wall  temperatures  roughly  coincide. 

These  factors  are  not  considered  too  serious  in  a  study  that  is 
mainly  concerned  with  fuel  type,  because  they  apply  with  equal  force  to 
all  fuels.  The  fact  that  the  measured  and  calculated  values  of  T 

W 

follow  the  same  trend,  as  evidenced  by  Figs.  44  to  47,  tends  to  support 
the  validity  of  using  the  luminosity  factor  concept  as  a  convenient  means 
for  incorporating  fuel  hydrogen  content  into  the  'standard*  equation  for 
flame  emissivity.  Thus  Eq.  (41)  may  be  rewritten  as 

eg  =  1  -  exp 


-97.44  P.(%H,r2  (q 


*b)0‘5  V1'5 


(50) 
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SECTION  IX 


EMISSIONS 

It  is  widely  recognized  that  combustion-generated  pollution  is  a 
threat  to  the  environment,  and  regulative  standards  have  been  imposed 
to  limit  the  pollutant  emissions  discharged  by  aircraft  engines 
operating  within  or  near  airports  [25,  26]. 

Smoke  is  the  most  obvious  pollutant  from  gas  turbine  engines 
because  it  can  be  seen  with  the  naked  eye.  Other  pollutants  of 
importance  are  carbon  monoxide  (CO),  unburned  hydrocarbons  (UHC) 
and  the  oxides  of  nitrogen  (N0X) . 

The  concentration  levels  of  most  pollutants  of  interest  in 
gas  turbine  exhausts  can  be  related  directly  to  the  temperature,  time 
and  concentration  histories  that  exist  in  the  combustor.  These 
histories  vary  from  one  combustor  to  another  and,  for  any  given 
combustor,  with  change  in  operating  conditions.  The  nature  of  pollutant 
formation  is  such  that  the  concentrations  of  carbon  monoxide  and  unburned 
hydrocarbons  are  highest  6c  low  power  conditions  and  diminish  with 
increase  in  power.  In  contrast,  oxides  of  nitrogen  and  smoke  are  fairly 
insignificant  at  lower  power  settings  and  attain  maximum  values  at  the 
highest  power  conditions. 

1 .  Carbon  Monoxide 

If  the  primary  zone  of  a  gas  turbine  combustor  is  designed  to 
operate  fuel  rich,  then  large  amounts  of  CO  will  be  formed  due  to  lack  of 
the  oxygen  needed  to  complete  the  reaction  to  C02.  If,  however,  the 
primary  zone  mixture  strength  is  stoichiometric  or  moderately  fuel  lean, 
then  significant  amounts  of  CO  will  be  present  due  to  the  dissociation 
of  C02.  In  principle  it  should  be  possible  to  reduce  this  CO  to  a 
negligible  level  by  the  staged  admission  of  additional  air  downstream 
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of  the  primary  zone  to  achieve  a  gradual  reduction  in  burned  gas 
temperature. 

In  practice,  CO  emissions  are  found  to  be  much  higher  than  predicted 
from  equilibrium  calculations,  and  to  be  highest  at  low  power  conditions, 
when  peak  temperatures  are  relatively  low.  This  is  in  conflict  with  the 
predictions  of  equilibrium  theory,  and  suggests  that  much  of  the  CO  arises 
from  incomplete  combustion  of  the  fuel.  This  may  be  caused  by  one  or 
more  of  the  following: 

1.  Inadequate  burning  rates  in  the  primary  zone  due  to  too  low  a  fuel/ 
air  ratio  and/or  insufficient  residence  time. 

2.  Inadequate  mixing  of  fuel  and  air,  which  produces  some  regions  in 
which  the  m’xcure  strength  is  too  weak  to  support  combustion,  and 
othey">  in  which  over-rich  combustion  yields  high  local  concentrations 
of  C  . 

3.  Quenching  of  the  post  flame  products  by  entrainment  with  the  liner 
wall-cooling  air.  (Annular  chambers,  because  of  their  lower  surface/ 
volume  ratio,  generally  give  lower  CO  emissions  than  tubular  systems.) 

CO  is  relatively  resistant  to  oxidation,  and  in  many  practical 
systems  its  oxidation  is  'rate  detemhr • no'  with  respect  to  the  minimum 
residence  time  and  mean  temperature  needed  for  complete  combustion.  At 
high  temperatures  the  major  reaction  removing  CO  is 

CO  +  OH  ^  C09  +  H9 
*-  2  2 

This  is  a  fast  reaction  over  z  breed  temperature  range.  At  lower  tempera¬ 
tures  the  reaction 

CO  +  h2o  l  co2  +  h2 

is  important  as  a  means  of  removing  CO. 
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2.  Unburned  Hydrocarbons 


Unburned  hydrocarbons  include  fuel  which  emerges  at  the  combustor 
exit  in  the  form  of  droplets  or  vapor,  and  also  the  products  of  thermal 
degradation  of  the  parent  fuel  into  species  of  lower  molecular  weight, 
such  as  methane  and  acetylene.  They  are  normally  associated  with  poor 
atomization,  inadequate  burning  rates,  the  chilling  effects  of  film-cooling 
air,  or  any  combination  of  these.  Increase  in  engine  power  setting  reduces 
the  emission  of  unburned  hydrocarbons,  partly  by  improved  fuel  atomization, 
but  mainly  through  the  effects  of  higher  inlet  air  pressure  and  temperature 
which  together  enhance  chemical  reaction  rates  in  the  primary  combustion 
zone. 

In  general,  the  emissions  of  UHC  parallel  those  of  CO.  Any  factor 
that  raises  the  level  of  CO  emissions  will  tend  to  raise  UHC  emissions 
also.  Conversely,  any  combustor/nozzle  modifications  carried  out  for 
the  reduction  of  CO  will  usually  reduce  UHC  also.  For  this  reason  no 
further  consideration  will  be  given  to  unburned  hydrocarbons. 

3.  Oxides  of  Nitrogen 

Oxides  of  nitrogen,  of  which  the  predominant  compound  at  high 
emission  levels  is  nitric  oxide,  are  produced  by  the  oxidation  of 
atmospheric  nitrogen  in  high-temperature  regions  o*  the  flame.  The 
process  is  endothermic  and  proceeds  at  a  significant  rate  only  at 
temperatures  above  around  1800  K.  Thus,  in  contrast  to  CO  and  UHC,  NO 
arises  only  in  the  hot  central  regions  of  the  combustor,  and  levels  are 
highest  at  full  power  conditions. 

Nitric  oxide  can  be  produced  by  three  different  mechanisms: 

1.  "Thermal"  NO  -  by  oxidation  of  atmospheric  nitrogen  in  the  post¬ 
flame  gases 
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2.  "Prompt"  NO  -  by  high  speed  reactions  at  the  flame  front 

3.  "Fuel"  NO  -  by  oxidation  of  nitrogen  contained  in  the  fuel. 

Thermal  NO.  It  has  been  fairly  well  established  that  NO  formation 

in  combustion  processes  proceeds  by  the  Zeldovich  chain  mechanism 

0o  *  20 

0  +  N,  *  NO  +  N 
Z  ■*- 

N  +  O0  l  NO  +  0 

The  main  fuel -air  reaction  proceeds  quickly  in  the  lean  case,  merely  adding 
heat  to  the  mixture  and  playing  a  minor  role  in  NO  formation.  Free  oxygen 
from  the  equilibrium  dissociation  of  unburned  oxygen  molecules  initiates 
the  chain.  A  nitrogen  atom  is  liberated  which  reacts  with  oxygen  to 
form  NO.  Equilibrium  dissociation  of  nitrogen  molecules  is  not  achieved 
at  the  temperatures  encountered  in  a  gas  turbine  combustor,  and  the  only  source 
of  nitrogen  atoms  is  the  second  reaction.  The  calculated  equilibrium  NO 
level  rises  continually  as  equivalence  ratio  is  reduced  at  a  fixed 
temperature,  or  as  temperature  is  increased  at  a  fixed  equivalence  ratio. 

In  a  practical  system, temperature  is  not  independent  of  equivalence  ratio, 
and  NO  formation  is  then  found  to  peak  on  the  fuel-lean  side  of 
stoichiometric.  This  is  a  consequence  of  competition  between  fuel  and 
nitrogen  for  the  available  oxygen.  Although  system  temperature  is  a 
maximum  at,  or  slightly  on  the  rich  side  of  stoichiometric,  the  available 
oxygen  is  then  all  consumed  preferentially  by  the  fuel  (due  to  the  higher 
rate  of  the  exothermic  fuel -oxygen  reaction).  At  equivalence  ratios 
below  around  0.8  the  fall-off  in  temperature  is  so  great  as  to  override 
the  effect  of  increasing  free  oxygen  concentration,  and  NO  levels  begin 
to  fall. 
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The  strong  correlation  between  temperature  and  N0X  emissions  was 
first  demonstrated  by  Lipfert  [27]  for  a  large  number  of  engines. 

Lipfert  based  his  correlation  on  combustor  inlet  temperature,  which  is 
important  only  insofar  as  it  influences  the  flame  temperature.  Roffe  and 
Venkataramani  [28]  have  examined  the  effects  of  combustion  temperature 
and  pressure  on  N0X  emissions,  using  propane/air  mixtures.  Their  results 
conform  well  with  the  expression 

ln(N0EI/t)  «  -72.28  +  2.80  ST  -  T/38  (51) 

where  T  is  the  adiabatic  flame  temperature,  K,  and  t  is  the  combustor 
residence  time,  ms.  They  found  no  effect  of  pressure  over  the  range  from 
0.5  to  3  MPa.  This  is  in  accordance  with  expectations,  since  the 
reaction  N2  +  02  =  2N0  involves  no  volume  change. 

The  effect  of  residence  time  on  N0X  emissions  was  examined  by 
Anderson  [29]  using  a  premixed/prevaporized  combustor  supplied  with 
gaseous  propane  fuel.  His  results  show  that  N0X  emissions  increase  with 
increase  in  residence  time  except  for  very  lean  mixtures  (<t>  =  0.4)  where 
the  rate  cf  formation  is  so  low  that  it  becomes  fairly  insensitive  to  time. 

Prompt  NO.  Under  certain  conditions,  and  especially  in  low  temperature, 
fuel-rich  flames,  NO  is  found  very  early  in  the  flame  region,  which  is 
in  conflict  with  the  idea  of  a  kinetically  slow  overall  process.  The 
mechanisms  involved  are  not  yet  fully  understood,  but  it  has  been  shown 
that  the  enhanced  reaction  rates  are  a  result  of  interactions  between 
the  many  intermediate  species  that  are  produced  during  the  main  hydro¬ 
carbon  -  CO  -  C02  reactions  [30]. 

Although  prompt  NO  levels  cannot  be  predicted  with  any  degree  of 
precision,  for  the  conditions  that  exist  in  modern  combustors  it  is  likely 


that  they  will  be  between  0  to  30  ppm  V,  with  low  values  occurring  at  high 
temperature,  lean  combustion  and  the  higher  values  at  low  temperature,  fuel 
rich  conditions. 

Fuel  NO.  If  fuels  contain  organically  bounded  nitrogen,  then  some  of 
this  nitrogen  will  eventually  form  the  so-called  "fuel  MO."  The  percentage 
of  nitrogen  undergoing  this  change  depends  on  the  nature  of  the  combustion 
process.  With  light  distillate  fuels,  the  content  of  organic  nitrogen  is 
small,  usually  less  than  0.06  percent.  Only  with  heavy  distillates  is 
the  content  of  organic  nitrogen  high  enough  to  yield  significant  amounts 
of  fuel  NO. 

4.  (J02  Emissions 

The  NO  produced  in  combustion  is  oxidized  to  NO2  as  soon  as  the  low 
temperatures  required  for  this  reaction  are  reached  after  leaving  the 
engine.  The  transformation  from  NO  into  N02  actually  begins  in  the 
combustor  in  zones  of  considerable  excess  air.  At  full  load  the  fraction 
of  N02  formed  within  the  combustor  is  very  small  but,  at  idle  conditions 
the  N02  level  may  be  as  much  as  50  percent  of  the  total  NOx  (NO  +  N02) 
emissions  [31,  32]. 

5.  Smoke 

Exhaust  smoke  is  caused  by  the  production  of  finely-divided  soot 
particles  in  fuel-rich  regions  of  the  flame,  and  may  be  generated  in  any 
part  of  the  combustion  zone  where  mixing  is  inadequate.  With  pressure 
atomizers  the  main  soot- forming  region  lies  inside  the  fuel  spray  at  the 
center  of  the  combustor.  This  is  the  region  in  which  the  recirculating 
burned  products  move  upstream  toward  the  fuel  spray,  and  where  local 
pockets  of  fuel  vapor  are  enveloped  in  oxygen-deficient  gases  at  high 
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temperature.  In  these  fuel-rich  regions  soot  may  be  produced  in  consider¬ 
able  quantities. 

Most  of  the  soot  produced  in  the  primary  zone  is  consumed  in  the  high 
temperature  regions  downstream.  Thus  from  a  smoke  viewpoint  a  cumbustor 
may  be  considered  to  comprise  two  separate  zones,  first,  the  primary  zone, 
which  governs  the  rate  of  soot  formation  and,  second,  the  intermediate 
zone  (and,  on  modern  high  temperature  engines,  the  dilution  zone  also) 
which  determines  the  rate  of  soot  consumption.  The  soot  concentration 
actually  observed  in  the  exhaust  gases  is  an  indication  of  the  dominance 
of  one  zone  over  the  other. 

Soot  is  not  an  equilibrium  product  of  combustion  except  at 
mixture  strengths  far  richer  than  those  employed  in  the  primary  zones  of 
combustors.  Thus  it  is  impossible  to  predict  its  rate  of  formation  and 
final  concentration  from  kinetic  or  thermodynamic  data.  In  practice, the 
rate  of  soot  formation  tends  to  be  governed  more  by  the  physical  processes 
of  atomization  and  fuel-air  mixing  than  by  kinetics. 

Many  specific  mechanistic  models  for  soot  formation  have  been  proposed. 
Current  thinking  tends  to  favor  the  notion  that  condensed  ring  aromatic 
hydrocarbons  may  produce  soot  via  a  different  mechanism  than  do  aliphatic 
hydrocarbons.  Aromatic  hydrocarbons  can  produce  soot  via  two  mechanisms: 

(1)  condensation  of  the  aromatic  rings  into  a  graphite-like  structure,  or 

(2)  breakup  to  small  hydrocarbon  fragments  which  then  polymerize  to  form 
larger,  hydrogen-deficient  molecules  which  eventually  nucleate  and  produce 
soot.  Based  on  his  shock  tube  studies,  Graham  [33,  34]  concludes  that  the 
condensation  route  is  much  faster  than  the  fragmentation/polymerization 
route.  According  to  this  simple  model,  aliphatics  produce  soot  via  the 
fragmentaticn/polymerization  mechanism  only:  As  a  result,  these  hydrocarbons 
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do  not  font  the  quantities  of  soot  produced  by  the  aromatics.  Indeed, 
during  the  fuel -rich  combustion  of  a  fuel  blend  composed  of  aromatics  and 
aliphatics,  the  aromatic  hydrocarbons  would  produce  the  major  quantity 
of  soot.  Combustion  of  the  aliphatic  portions  of  the  fuel  would  influence 
temperature  and  hydrocarbon  fragment  concentration  but  soot  formation  via 
fragmentation/polymerization  would  be  minimal. 

Experimental  data  obtained  by  Blazowski  [35]  using  various  blends 
of  iso-octane  and  toluene  fuels  were  found  to  be  consistent  with  this 
model.  However,  the  results  of  an  experimental  study  by  Naegeli  and 
Moses  [36]  suggest  that  the  picture  will  be  more  complicated  for  fuels 
with  high  concentrations  of  polycyclic  aromatics. 

For  gas  turbine  combustors  the  main  controlling  factors  for  soot 
formation  and  smoke  have  been  determined  experimentally  as  fuel  properties, 
combustion  pressure  and  temperature,  fuel/air  ratio,  atomization  quality 
and  mode  of  fuel  injection  [37  to  50]. 

Influence  of  Fuel  Properties.  Fuel  properties  can  influence  smoke 
production  in  two  ways:  firstly,  by  inducing  the  formation  of  local  over¬ 
rich  fuel  regions  and,  secondly,  by  exerting  variable  resistance  to  carbon 
formation.  The  former  is  controlled  by  physical  properties  such  as 
viscosity  and  volatility  which  affect  the  mean  drop  size,  penetration  and 
rate  of  evaporation  of  the  fuel  spray,  whereas  the  latter  relates  to 
molecular  structure.  It  is  well  established  that  smoking  tendency 
increases  with  reduction  in  hydrogen  content  [51].  However,  it  is  of 
interest  to  note  that  Naegeli  and  Moses  [36]  found  that  fuels  containing 
high  concentrations  of  polycyclic  aromatics  produced  more  soot  than  would 
be  expected  from  a  correlation  based  solely  on  hydrogen  content. 


Influence  of  Pressure.  Problems  of  soot  and  smoke  are  always  most 


severe  at  high  pressures.  There  are  several  reasons  for  this,  some  of 
which  derive  from  chemical  effects,  while  others  stem  from  physical  factors 
which  affect  spray  characteristics  and  hence  also  the  distribution  of 
mixture  strength  in  the  soot-generating  regions  of  the  flame  [38,  43-45]. 

For  premixed  kerosine-air  flames  it  is  found  that  no  soot  is  formed  at 
pressures  below  around  0.6  MPa,  or  at  equivalence  ratios  below  about  1,3  [42]. 
One  adverse  effect  of  increase  in  pressure  is  to  extend  the  limits  of 
flammability,  so  that  soot  is  produced  in  regions  which,  at  lower  pressures, 
would  be  too  rich  to  burn.  Increase  in  pressure  also  accelerates 
chemical  reaction  rates,  so  that  combustion  is  initiated  earlier  and  a 
larger  proportion  of  the  fuel  is  burned  in  the  fuel-rich  regions  adjacent 
to  the  spray.  With  pressure  atomizers  of  the  duplex  or  dual-orifice  type, 
reduced  spray  penetration  is  one  of  the  main  causes  of  smoke  at  high 
pressures  [45].  Thus  whereas  at  low  pressures  the  fuel  is  distributed 
across  the  entire  combustion  zone,  at  high  pressures  it  tends  to  concentrate 
in  the  soot-forming  region  just  downstream  of  the  nozzle.  Another  adverse 
effect  of  an  increase  in  pressure  is  to  reduce  the  cone  angle  of  the  spray. 
This  encourages  soot  formation,  partly  by  increasing  the  mean  drop  size, 
but  mainly  by  raising  the  mixture  strength  in  the  soot-forming  zone.  The 
total  effect  of  all  these  factors  is  that  smoke  emission  increases  steeply 
with  pressure.  Airblast  atomizers  are  spared  these  problems  because  the 
fuel  drops  they  produce  are  always  airborne  and  their  distribution  throughout 
the  combustion  zone  is  dictated  solely  by  the  liner  airflow  pattern  which 
is  not  susceptible  to  changes  in  pressure. 


Influence  of  Temperature.  Combustor  inlet  temperature  tends  to 
influence  many  variables  in  the  carbon  formation  process,  so  that  the 
effect  of  a  change  in  on  exhaust  smoke  is  not  always  clear. 

However,  it  is  well  established  that  increase  in  combustor  outlet 
temperature  reduces  smoke  by  extending  the  soot  consumption  region  down¬ 
stream  into  the  dilution  zone  where  there  is  ample  air  to  bum  up  the 
soot,  provided  there  is  sufficient  temperature  also. 

Influence  of  Fuel/Air  Preparation.  Soot  is  formed  only  in  fuel -rich 
regions  of  the  flame.  Thus  all  that  is  needed  to  eliminate  soot  and  smoke 
is  to  ensure  that  nowhere  in  the  flame  region  does  the  equivalence  ratio 
exceed  around  1.3.  However,  even  when  the  overall  equivalence  ratio  in 
the  primary  zone  is  well  below  1.3,  imperfections  in  fuel-air  mixing  can 
create  local  regions  in  which  pockets  of  fuel-rich  mixture  are  enveloped 
in  oxygen-deficient  gases  at  high  temperature,  leading  to  high  rates  of 
soot  formation.  The  superior  performance  of  the  airblast  atomizer  is  due 
partly  to  better  atomization  but  also  to  the  thorough  fuel-air  mixing 
incurred  in  the  atomization  process  prior  to  combustion,  which  effectively 
eliminates  fuel-rich  pockets  from  the  combustion  zone. 

In  practice  it  is  found  that,  even  at  the  worst  conditions,  only  a 
very  small  fraction  of  the  fuel  carbon  is  converted  to  soot,  and  almost 
all  of  this  is  oxidized  in  the  regions  downstream  of  the  primary  combust-'  r: 
zone  [44].  Appleton  [52]  has  developed  a  method  for  calculating  the 
size  reduction  experienced  by  a  carbon  particle  as  it  travels  through  a 
mixture  of  known  pressure,  temperature  and  composition.  For  conditions 
typical  of  the  gas  turbine,  this  method  predicts  that  all  particles  below 
0.04  um  will  be  oxidized,  and  that  the  maximum  rate  of  oxidation  proceeds 
at  an  equivalence  ratio  of  approximately  0.75.  A  similar  conclusion  was 
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reached  by  Fenimore  and  Jones  [53]  who  found  evidence  that  as  soot 
aggregates  increase  in  size  they  become  more  resistant  to  oxidation. 

According  to  these  workers  a  soot-bearing  flame  ought  to  burn  smoke-free 
if  the  carbon  aggregates  are  "less  than  0.01  ym,  no  matter  how  many  are 
present. 

6.  Prediction  of  Emissions  Characteristics 

Many  attempts  have  been  made  to  model  gas  turbine  combustors  with  a 
view  to  predicting  their  emission  characteristics.  Most  modeling  has  been 
concerned  with  NQV  emissions,  but  efforts  have  also  been  made  to  predict 

A 

the  formation  of  other  important  pollutant  species.  To  be  successful  a 
model  must  accommodate  the  complex  flow  behavior  and  include  a  kinetic 
scheme  of  the  important  chemical  reactions  occurring  within  the  combustor. 

The  kinetics  of  some  relevant  combustion  processes  are,  unfortunately, 
not  well  understood  at  the  present  time,  particularly  for  the  production 
of  carbon,  carbon  monoxide  and  the  hydrocarbon  species  that  are  intermediaries 
in  the  fuel  oxidation  process. 

The  primary  requirement  for  a  satisfactory  emissions  model  for  gas- 
turbine  combustors  is  that  it  should  represent  an  optimum  balance  between 
accuracy  of  representation,  utility,  ease  of  use,  economy  of  operation, 
and  capability  for  further  improvement.  In  recent  years  considerable 
effort  has  been  directed  toward  the  development  of  relatively  complex 
mathematical  emissions  models  that  can  be  applied  to  gas  turbines.  They 
have  varied  in  level  of  sophistication  from  those  with  potential  to  yield 
a  complete  description  of  the  relevant  thermodynamic  and  chemical  proper¬ 
ties  as  a  function  of  spatial  location  within  the  combustor,  to  others 
which  merely  assume  that  homogeneous  conditions  exist  at  all  axial 
stations  [54-63]. 
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The  high  cost  and  complexity  of  the  more  sophisticated  mathematical 
models,  combined  with  the  somewhat  dubious  value  of  the  results  obtained, 
have  encouraged  the  development  of  semi -empirical  models  for  N0X  and  CO 
emissions.  Hung's  approach  to  the  modeling  of  N0X  emissions  places  more 
emphasis  on  the  physical  processes  considered  tc  be  important  and  delegates 
chemical  kinetics  to  a  relatively  minor  role  [64].  His  analysis  includes 
a  five-region  combustor  internal  flowfield  model,  fuel  distribution 
models  for  liquid  and  gaseous  fuels,  a  single  overall  hydrocarbon  complete 
combustion  model,  a  nitric  oxide  formation  model  based  on  the  Zeldovich 
mechanism,  a  diffusion-limited  complete  mixing  model,  and  a  model  to 
account  for  the  influence  of  ambient  humidity.  This  model  has  been  used 
successfully  in  predicting  the  influence  on  N0X  emissions  of  water 
injection  and  wide  variations  in  fuel  type  [64-67]. 

Other  successful  semi-empirical  models  for  predicting  emissions  have 
been  developed  by  Fletcher  and  Heywood  [55,  68]  and  by  Hanmond  and  Mellor 
[69-71].  Useful  critical  evaluations  of  both  mathematical  and  semi- 
empirical  prediction  methods  have  been  made  by  Rubins  and  Marchionna  [72], 
Sullivan  and  Mas  [73],  and  Odgers  [74]. 

Empirical  models  can  also  play  an  important  role  in  the  design  and 
development  of  low  emission  combustors.  They  may  serve  to  reduce  the 
complex  problems  associated  with  emissions  to  forms  which  are  more  meaning¬ 
ful  and  tractable  to  the  combustion  engineer,  who  often  requires  only  an 
insight  and  a  quick  estimate  of  the  levels  attainable  with  the  design 
variables  at  his  disposal.  They  also  permit  acre  accurate  correlations  of 
emissions  for  any  one  specific  combustor  than  can  be  achieved  by  the  more 
general,  analytical  models  discussed  previously. 
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7.  N0X  and  CO 

In  attempting  to  derive  an  empirical  model  for  emissions,  emphasis  is 
placed  on  CO  and  N0x-  This  is  because  the  highly  complex  and  unknown 
nature  of  the  hydrocarbon  oxidation  reaction  makes  it  impossible  to  derive 
satisfactory  models  for  UHC  and  smoke.  Moreover,  it  is  known  that  reduction 
in  CO  and  N0x  will  also  ensure  lower  levels  of  UHC  and  smoke,  respectively. 

It  may  be  assumed  for  both  nitric  oxide  and  carbon  monoxide  that 
their  exhaust  concentration  is  proportional  to  the  product  of  three 
terras  which  are  selected  to  represent  the  following: 

1.  Mean  residence  time  in  the  combustion  zone. 

2.  Chemical  reaction  rates,  and 

3.  Mixing  rates. 


Expressions  for  these  three  parameters  may  be  derived  in  simplified  form 
as: 


1 . 


Residence  time  =  k  = 

mA  mART 
(  \ 


i.e.,  residence  time  « 


PV 

y 


2.  It  is  assumed  that  reaction  rates  are  a  function  of  pressure 
and  temperature  only,  i.e. 

reaction  rate  <*  Pmexp  (zT)  for  N0X 

and  reaction  rate  Pnexp  (cT)  for  CO. 


3.  It  is  assumed  that  mixing  rates  are  a  function  of  liner 
pressure  drop.  Specifically  we  have 


mixing  rate  = 


f&P)X 
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Thus  NQX  =  -j-  (residence  time)  (reaction  rate)  (mixing  rate) 
or 
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(52) 


where  y  =  1  +  m,  and  A  is  a  constant. 

Similarly,  for  CO  we  have 

CO  =  -f  (residence  time)-1 (reaction  rate) (mixing  rate) 


or 


CO  =  c 


f. 
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AP 

PV 
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cj 

,  } 

Pnexp(-cT) 


=  C  Vc_,mA  T  (AP/P)aP  bexp(-cT) 


(53) 


where  b  -  n  -  1,  and  C  is  a  constant. 

It  is  recognized  that  the  above  equations  have  no  strong  theoretical 
foundation.  However,  they  do  embody  the  main  variables  of  combustor  size, 
pressure  loss,  flow  proportions  and  operating  conditions  of  inlet  pressure, 
temperature  and  air  mass  flow.  The  effect  of  variations  in  overall 
combustor  fuel/air  ratio  is  also  included  via  its  influence  on  primary- 
zone  temperature.  Fuel  type  affects  both  flame  temperature  and  mean 
drop  size.  For  N0X  drop  size  is  unimportant  since,  at  the  high 
pressure  conditions  where  N0X  emissions  are  most  prominent,  the 
fraction  of  the  total  combustion  volume  employed  in  fuel  evaporation  is 
so  small  that  wide  variations  in  fuel  drop  s’ize  hove  a  negligible 
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effect  on  N0X  emissions.  However,  at  low  pressure  operation  where 
CO  emissions  attain  their  highest  concentrations,  a  significant  proportion 
of  the  primary-zone  volume  is  needed  to  evaporate  the  fuel.  Under  these 
conditions,  any  factor  that  influences  fuel  evaporation  rates,  such  as 
evaporation  constant,  or  mean  drop  size,  will  have  a  direct  effect 
on  the  volume  available  for  chemical  reaction  and,  therefore,  on  the 
emissions  of  CO  and  UHC.  Thus,  for  the  correlation  of  CO  data  the 
effects  of  fuel  type  cannot  be  ignored.  The  manner  in  which  they  may  be 
introduced  conveniently  into  the  equation  for  CO  emissions  is  illustrated 
below.  From  analysis  of  the  experimental  data  contained  in  references 
[1  to  6J  it  was  found  for  Eq.  (52)  that 

A  =  9  x  10“8,  x  =  0,  y  =  1.25,  z  =  0.01. 

For  Eq.  (53)  the  results  show  that 

C  -  H6,  a  =  -0.5,  b  *  -1.5,  c  =  0.00345 

Substituting  these  values  into  Eqs.  (52)  and  (53)  gives 


9  x  10"8  pJ-25  V  exp  (0.01  T  , ) 

N0  =  - L - 2 - g/kg 

mA  Tpz 


(54) 
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From  Eq.  (54)  it  may  be  noted  that  the  only  influence  of  fuel  type 
on  N0X  formation  is  via  the  two  temperature  terms  Tp?  and  Tsfc,  The 


former  is  calculated  as 
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pz 
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where  ATp2  is  the  temperature  rise  due  to  combustion  corresponding  to 
the  inlet  temperature,  Tg,  and  fuel /air  ratio,  (qQV/f).  T$t  is  the 
stoichiometric  flame  temperature  corresponding  to  the  inlet  temperature, 

Tg.  Equation  (54)  suggests  that  in  the  combustion  of  heterogeneous 
fuel -air  mixtures  it  is  the  stoichiometric  flame  temperature  that 
determines  the  formation  of  N0X-  However,  for  the  residence  time  in 
the  combustion  zone,  which  is  also  significant  to  N0X  formation,  the 
appropriate  temperature  term  is  the  bulk  value,  Tpz,  as  indicated  in 
the  denominator  of  Eq.  (54). 

Equation  (54)  is  suitable  for  conventional  spray  combustors  only. 

For  lean  premix/prevaporize  combustors,  in  which  the  maximum  attainable 
temperature  is  Tpz,  Eq.  (54)  may  still  be  used,  provided  that  Tpz  is 
substituted  for  T  It  should  also  be  noted  that  predictions  of  N0X 
based  on  Eq.  (54)  tend  to  be  too  high  when  the  overall  combustor  air/fuel 
ratio  exceeds  a  value  of  around  100.  This  is  because  with  diminishing 
fuel/air  ratio  the  flame  shrinks  back  toward  the  fuel  nozzle  and  no  longer 
occupies  the  entire  combustion  volume,  Vc-  However,  this  is  not  considered 
a  serious  drawback  since,  in  practice,  interest  is  normally  focused  on 
conditions  of  high  fuel/air  ratio,  where  NQX  formation  rates  attain  their 
highest  values. 

The  excellent  correlation  of  experimental  data  on  NO  provided  by 

A 

Eq.  (54)  is  illustrated  in  Figs.  48  to  58.  These  figures  include  all 
combustors  except  the  J85  for  which  the  measured  values  are  too  low  for 
satisfactory  correlation. 

The  formation  of  CO  in  the  primary  combustion  zone  takes  appreciably 
longer  than  the  time  required  to  produce  N0X-  In  consequence,  the 
relevant  temperature  is  not  the  local  peak  value  adjacent  to  the  evaporating 
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Figure  51.  Comparison  of  Measured  and  Predicted  Values  of  NO 
Emissions  for  J79-17C  Combustor.  (Fuels  7A  to  13ft) 
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gure  52,  Comparison  of  Measured  and  Predicted  Values  of  NO 
Emissions  for  F101  Combustor.  (Fuels  1  to  6). 
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Figure  54.  Comparison  of  Measured  and  Predicted  Values  of  NO 
Emissions  for  TF  41  Combustor.  (Fuels  I  to  6). 
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Figure  55.  Comparison  of  Measured  and  Predicted  Values  of  NO 
Emissions  for  TF  41  Combustor.  (Fuels  7  to  12).  x 
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Figure  56.  Comparison  of  Measured  and  Predicted  Values  of 
N0„  Emissions  for  TF39  Combustor. 
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Figure  57.  Comparison  of  Measured  and  Predicted  Values  of  NO 
Emissions  for  TF  33  Combustor.  x 
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Ail  Fuels 


fuel  drops,  but  the  average  value  throughout  the  primary  zone,  namely  7  . 

pz 

Also,  because  CO  emissions  are  most  important  at  low  pressure  conditions, 
where  evaporation  rates  are  relatively  slow,  it  is  necessary  to  reduce 
the  combustion  volume,  Vc,  by  the  volume  occupied  in  fuel  evaporation,  Vg. 
This  is  given  by  (see  Pattern  Factor  Section) 

*e  ■  °-55  fpz  "A  °o'V  ‘eff  <“> 


It  is  of  interest  to  note  the  inclusion  of  a  pressure  loss  term  in 
Eq.  (55).  This  suggests  that  the  higher  turbulence  created  by  an  increase 
in  liner  pressure  drop  promotes  better  mixing  in  the  combustion  zone 
and  helps  to  eliminate  the  rich  and  weak  pockets  of  fuel-air  mixture, 
both  of  which  are  conducive  to  high  rates  of  CO  formation. 

The  very  satisfactory  correlation  of  experimental  data  on  CO 
emissions  obtained  with  Eq.  (55)  is  illustrated  in  Figs.  59  to  63, 
for  J79-17A,  J79-17C,  F101,  TF41  and  F100  combustors  respectively. 


8.  Smoke 

In  order  to  express  r  noke  emissions  in  the  same  basic  units  as 
those  employed  for  N0X  and  CO,  the  first  step  must  be  to  convert  the 
smoke  numbers  (SN)  quoted  in  references  [1  to  6]  to  soot 
concentrations  (Xc)  expressed  in  mg/kg.  This  conversion  was  accomplished 
using  the  following  different  factors  for  different  levels  of  stroke 
number  [75]. 


SN  =  0  to  1 

SN  =  1  to  5 
SN  =  5  to  10 
SN  -  10  to  20 
SN  =  20  to  30 
SN  >  30 


Xc  =  0.1  (SN) 

log  Xc  =  0.136  (SN)  -  1.136 
log  Xc  =  0.06265  (SN)  -  0.769 
log  Xc  =  0.03187  (SN)  -  0.4614 
log  Xc  =  0.0301  (SN)  -  0.426 
log  Xc  =  0.02538  (SN)  -  0.2845 
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Figure  60.  Comparison  of  Measured  and  Predicted  Values  of  CO 
Emissions  for  J79-17C  Combustor. 
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Figure  61.  Comparison  of  Measured  and  Predicted  Values  of  CO 
Emissions  for  FI01  Combustor. 
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62.  Comparison  of  Measured  and  Predicted  Values  of  CO 
Emissions  for  TF  41  Combustor. 
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Figure  63.  Comparison  of  Measured  and  Predicted  Values  of  CO 
Emissions  for  FIOO  Combustor. 
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The  following  equation  was  used  to  convert  engine  exhaust  soot 
concentrations  into  corresponding  combustor  exit  values. 


X  =  X  Nl  +  q8 
c4  c8  q8  1  +  q4 


(57) 


For  the  purpose  of  analysis,  it  is  convenient  to  consider  two  separate 
zones  (a)  a  soot-forming  zone,  and  (b)  a  soot  oxidation  zone.  The  soot 
concentration  measured  at  the  combustor  exit  represents  the  difference  in 
effectiveness  between  these  two  competing  processes.  Unfortunately, 
any  attempt  to  derive  suitable  expressions  to  represent  rates  of  soot- 
formation  and  soot-oxidation  is  seriously  hampered  by  lack  of  knowledge 
of  the  basic  mechanisms  involved,  so  that,  in  practice,  there  is  little 
alternative  except  to  resort  to  an  empirical  approach.  Useful  guidance 
is  provided  by  the  knowledge  gained  from  past  experience  in  attempting 
to  alleviate  the  problems  of  smoke  and  soot  formation  in  gas  turbine 
combustors.  Thus,  for  example,  the  work  of  Macfarlane  [42]  has 
shown  that  soot  formation  increases  rapidly  with  increase  in  pressure, 
and  is  appreciably  diminished  by  increase  in  AFR.  Moreover,  sufficient 
is  known  to  indicate  that  soot  oxidation  proceeds  most  rapidly  in 
regions  of  high  temperature  containing  excess  air.  These  considerations, 
in  conjunction  with  analysis  of  the  experimental  data,  lead  to  the 
following  expressions  for  the  suot  formation  and  soot  oxidation 
processes . 
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(59) 


[C  q  P,  exp  (0.0014  T  )]  (Aromatic  Content  +  B) 
_  P_  ^  S  i- 

*pz  ™A  ^sz  Tpz 


Xc  "  XF  ‘  Xo 


fpz  mA  Tpz 


A  -  --  exp  (0.0014  T  ) 
qsz 


Aromatic  Content  +  B 


(60) 


Application  of  this  equation  to  the  correlation  of  experimental  data  on 
soot  concentrations  contained  in  references  1  to  4  yields  the  results  shown 
in  Figs.  64  to  67.  The  values  of  A,  B  and  C  associated  with  these  figures 
are  listed  in  Table  6. 

If  allowance  is  made  for  the  difficulties  involved  in  the  sampling 
and  measurement  of  soot  concentrations,  and  the  poor  measurability  of  fuel 
aromatics  content,  the  level  of  agreement  between  measured  and  predicted 
values  of  soot  concentration,  as  exhibited  in  Figs.  S4  to  67,  is  clearly 
very  reasonable.  However,  if  these  results  are  accepted  at  face  value  they 
could  also  be  very  misleading.  Although  Eq.  (60)  predicts  quite  well  the 
influence  of  combustor  operating  conditions  on  smoke  output,  and  also 
demonstrates  that  soot  concentrations  rise  with  increase  in  aromatic  content 
of  the  fuel,  it  also  shows  that  the  extent  of  this  increase  varies  from 
one  combustor  to  another  in  a  manner  that  cannot  be  predicted  a  priori. 

Thus  Eq.  (60)  offers  no  guidance  whatsoever  on  the  likely  smoke  emissions 
to  be  anticipated  from  any  new  type  of  combustor.  Only  if  the  values  of  A, 

B  and  C  remained  sensibly  constant  for  all  combustors  would  it  be  reasonable 
to  regard  Eq.  (60)  as  completely  satisfactory  for  the  prediction  of  smoke 
emissions.  In  this  context  it  is  of  interest  to  note  that  the  values  of 
A  and  C  for  the  079-17C  and  rF41  combustors  are  almost  identical.  However, 
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Figure  64.  Graphs  Illustrating  Influence  of  Aromatics  Content 

and  E^ne  Operating  Conditions  on  Soot  Emissions 
for  079- 1 7A  Combustor. 
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Figure  66.  Graphs  Illustrating  Influence  of  Aromatics  Content 
and  Engine  Operating  Conditions  on  Soot  Emissions 
for  FI 01  Combustor. 


Figure  67.  Graphs  Illustrating  Influence  of  Aromatics  Content 

and  Engine  Operating  Conditions  on  Soot  Emissions 
for  TF41  Combustor. 


TABLE  6.  Values  of  A,  B  and  C  employed  in  Eq.  (60). 
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ENGINE 

A 

B 

C 

J79-17A 

32.9 

-0.06 

0.08 

J79-17C 

3.22 

-0.02 

0.0078 

F101 

0.72 

+0.55 

0.001 

TF41 

3.2 

+0.16 

0.007 
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for  the  other  two  combustors  examined  the  corresponding  values  are  widely 
different,  indicating  that  Eq.  (60)  fails  to  take  into  account  one  or  more 
processes  of  great  importance  to  soot  formation.  One  obvious  omission  is 
a  term  to  describe  the  degree  of  mixing  of  fuel  and  air  prior  to  combustion. 
This  is  known  to  have  a  strong  influence  on  smoke  output;  for  example,  the 
very  large  difference  in  smoke  output  between  the  J79-17A  and  J79-17C 
combustors,  as  reflected  in  the  ten-fold  difference  in  the  value  of  A,  is 
known  to  be  due  in  large  measure  to  the  steps  taken  to  improve  the  premixing 
of  fuel  and  air  in  the  latter  case.  Thus  Eq.  (60)  contains  some  serious 
deficiencies,  but  little  improvement  in  the  prediction  of  smoke  emission 
can  be  expected  until  more  quantitative  information  is  available  on  the 
influences  of  fuel -air  preparation  and  fuel  chemistry  on  soot  formation. 
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PATTERN  FACTOR 

Perhaps  the  most  important  and,  at  the  same  time,  most  difficult 
problem  in  the  design  and  development  of  gas  turbine  combustion  chambers, 
is  that  of  achieving  a  satisfactory  and  consistent  distribution  of  temper¬ 
ature  in  the  efflux  gases  discharging  into  the  turbine.  In  the  past, 
experience  has  played  a  major  role  in  the  determination  of  dilution-zone 
geometry,  and  trial  and  error  methods  have  been  widely  used  in  developing 
the  temperature-traverse  quality  of  individual  combustor  designs  to  a 
satisfactory  standard.  Experimental  investigations  into  dilution-zone 
performance  carried  out  on  actual  chambers  have  provided  useful  guidance, 
but  very  often  it  has  proved  difficult  or  impossible  to  distinguish  the 
separate  influences  of  all  the  variables  involved.  Thus  although  it  is 
now  generally  accepted  that  a  satisfactory  temperature  profile  is  dependent 
upon  adequate  penetration  of  the  dilution  jets,  coupled  with  the  correct 
number  of  jets  to  form  sufficient  localized  mixing  regions,  the  manner  in 
which  the  total  dilution-hole  area  is  utilized  in  terms  of  number  and  size 
of  holes  is  still  largely  a  matter  of  experience.  Unfortunately,  more 
basic  studies  of  jet  mixing  do  not  usually  yield  results  that  can  readily 
be  expressed  in  terms  of  the  parameters  most  familiar  to  those  concerned 
with  combustion-chamber  design.  However,  some  of  these  investigations  can 
provide  a  useful  guide  to  the  relationships  involved. 

The  temperature  obtained  by  an  elemental  volume  of  gas  at  the  chamber 
outlet  is  dependent  on  the  history  It  undergoes  from  the  time  It  emerges 
from  the  compressor.  During  Its  passage  through  the  combustor  Its  temper¬ 
ature  and  composition  changes  rapidly  under  the  Influence  of  various 
combustion,  heat  transfer  and  mixing  processes,  none  of  which  are  perfectly 
understood.  Considering  only  the  final  mixing  process,  this  is  affected 
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in  a  complicated  manner  by  the  dimensions,  geometry,  and  pressure  drop  of 
the  liner,  the  size,  shape  and  discharge  coefficients  of  the  liner  holes, 
the  airflow  distribution  to  various  zones  of  the  chamber,  and  the  temperature 
distribution  of  the  hot  gases  entering  the  dilution  zone.  For  any  given 
combustor,  the  latter  is  strongly  influenced  by  fuel  spray  characteristics 
such  as  drop  size,  spray  angle  and  spray  penetration,  since  these  control 
the  pattern  of  burning  and  hence  the  distribution  of  temperature  in  the 
primary-zone  effux.  It  is  known  that  spray  characteristics  are  strongly 
influenced  by  pressure,  especially  with  atomizers  of  the  simplex  or  dual¬ 
orifice  type,  and  it  is  to  be  expected  therefore  that  temperature  traverse 
will  also  vary  with  pressure,  although  the  extent  of  this  variation  will 
vary  from  one  chamber  to  another  depending  on  design  and,  in  particular, 
on  length. 

The  most  important  temperature  parameters  are  those  that  affect  the 
power  output  of  the  engine  and  the  life  and  durability  of  the  hot  sections 
downstream.  As  far  as  the  overall  engine  performance  is  concerned  the  most 
important  temperature  is  the  turbine  inlet  temperature.  ,  which  is  the 
mass  flow  weighted  mean  of  the  combustor  exit  t.e  erature.  Since  the 
nozzle  guide  vanes  are  fixed  relative  to  the  combustor  they  must  be 
designed  to  withstand  the  maximum  temperature  found  in  the  traverse. 

Thus,  the  parameter  of  most  relevance  to  the  design  of  nozzle  guide  vanes 
is  the  overall  temperature  distribution  factor,  which  highlights  this 
maximum  temperature.  It  is  normally  defined  as:  - 

T  -  T. 

Pattern  factor  =  -2^ — __  (51) 

‘  3 
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The  temperatures  of  most  significance  to  the  turbine  blades  are  those 
that  constitute  the  average  radial  profile.  This  is  obtained  by  adding 
together  the  temperature  measurements  around  each  radius  of  the  liner  and 
then  dividing  by  the  number  of  locations  at  each  radius*  i.e.,  calculating 
the  arithmetic  mean  at  each  radius.  The  expression  used  to  describe  the 
radial  temperature  distribution  factor,  also  known  as  the  profile  factor, 
is 

Profile  factor  =  -5^ — -i-  (62) 

4  "  '3 


In  Eqs.  (61)  and  (62) 

T|nax  =  maximum  recorded  temperature 

T  =  maximum  circumferential  mean  temperature 
mr 

T^  =  mean  inlet  air  temperature 
T4  =  mean  exit  temperature 
1 .  Correlation  of  Data 

Two  parameters  of  crucial  importance  to  pattern  factor  are  liner 
length,  which  controls  the  time  and  distance  that  are  available  for  mixing, 
and  the  pressure  drop  across  the  liner,  which  governs  the  penetration  of 
the  dilution  jets  and  their  rate  of  mixing  with  the  products  of  combustion. 
From  analysis  of  experimental  data  on  tubular,  tubo-annular  and  annular 
combustors,  it  is  found  that  [9] 


/ 


V 

qref 


(63) 
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where  Af>|/qref  =  ^iner  Pressure  loss  factor 
=  total  liner  length 
Dl  =  liner  diameter,  or  height 


The  data  correlations  obtained  for  tubular  and  annular  liners  are  shown 
in  Figs.  68  and  69  respectively.  In  connection  with  these  figures  it  should 
be  noted  that  the  correlation  is  based  not  on  the  L/D  ratio  of  the  dilution 
zone,  but  on  that  of  the  complete  liner.  This  is  found  to  provide  a  better 
fit  to  the  data. 

For  tubular  and  tubo-annular  combustors  it  is  found  that 


<Tmax  '  -  T^>  "  1  -  e*P 


0.070  (LL/DL)(APL/qref) 


-1 


(64) 


while  for  annular  combustors 


(Tmax  ‘  ty'<T4  -  T3>  =  1  *  e*P  ‘  °’050 


-1 


(65) 


Comparison  of  Figs.  68  and  69  reveals  that,  for  any  given  value  of 
(LL/DL) ,  the  pattern  factor  of  the  tubular  system  is  superior  to  that  of  an 
annular  system.  This  is  because  tubo-annular  combustors  tend  to  have  lower 
annulus  velocities  than  annular  combustors,  a  property  that  is  conducive  to 
better  jet  penetration  and  mixing.  Moreover,  with  tubo-annular  combustors 
the  pressure  distributions  in  the  annulus  are  less  susceptible  to  changes 
in  inlet  velocity  profile.  This  is  because  the  longitudinal  gap  between 
liners  ensures  automatic  radial  balancing  of  pressure  along  the  entire 
combustor  length.  In  marked  contrast,  in  the  annular  combustor  the 
pressure  difference  between  the  inner  and  outer  annuli  is  determined  by 
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Figure  68.  Pattern  Factor  Correlation  for  Tube- Annular  Combustors. 


Figure  69.  Pattern  Factor-  Correlation  for  Annular  Combustors 


the  shape  of  the  total  pressure  profile  in  the  flow  when  it  reaches  the 
snout. 

2.  Data  Analysis 

As  no  information  on  liner  pressure  loss  factor  is  contained  in 
references  1  to  6  this  might  appear  to  rule  out  the  use  of  Eqs.  (64)  and 
(65)  in  the  present  analysis.  Another  drawback  to  these  equations  is 
the  inherent  assumption  that  all  the  liner  length  is  utilized  in  mixing 
and  combustion  and  that  the  length  occupied  by  evaporation  processes  is 
essentially  zero.  Although  this  assumption  is  not  unreasonable  for 
highly  volatile  fuels  of  low  viscosity,  such  as  JP4,  it  is  difficult  to 
justify  for  some  of  the  alternative  fuels  employed  in  this  program.  These 
deficiencies  may  be  remedied  to  a  large  extent  by  rewriting  Eqs,  (64) 
and  (65)  in  following  form 


exp  -  Q 


1 

AP, 

L 

[qref 

Ll-  Le 

,  7 

(66) 


where  Q  =  0.070  for  tubular  liners 
or,  Q  =  0.050  for  annular  liners 


Lg  is  the  liner  length  required  to  evaporate  the  fuel  spray.  This 
is  obtained  as  the  product  of  the  average  predilution  liner  velocity  and 
the  evaporation  time. 
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The  constant  of  0.6  in  Eq.  (67)  stems  from  the  assumption  that  60  percent 
of  the  total  combustor  airflow  enters  the  liner  upstream  of  the  liner.  A 
more  accurate  assessment  of  this  constant  for  each  individual  liner  design 
would  not  be  justified  at  this  stage  due  to  the  larger  uncertainty  sur¬ 
rounding  the  value  of  0Q.  In  practice  it  was  found  that  correlation  of 
experimental  data  could  be  improved  by  reducing  the  value  of  this  constant 
from  0.6  to  0.33.  This  may  be  justified  on  the  grounds  that  the  combustion 
process  does  not  wait  until  the  fuel  has  fully  evaporated;  instead  burning 
commences  as  soon  as  sufficient  fuel  has  evaporated  to  produce  a  flammable 
mixture.  Thus  Eq.  (66)  may  be  rewritten  as 


1  -  exp  -  Q' 


•  2)1  -1 
°-33  mA  V 

pg  AlAAeff 


(68) 


where  pg  is  the  average  gas  density  upstream  of  the  dilution  zone.  It  is 
calculated  at  a  temperature  Tg  which  is  obtained  as 


where  £>Tg  is  the  temperature  rise  due  to  combustion  for  a  fuel/air  ratio  of 
0.6  qQv.  Al  is  the  average  cross-sectional  area  of  the  liner.  It  is 
estimated  by  dividing  the  volume  of  the  liner  by  its  maximum  length. 
is  the  average  diameter  or  height  of  the  liner.  For  a  tubular  liner  it  is 
readily  obtained  as  =  (4  A^/ir)0'5. 

Mean  drop  sizes  for  insertion  into  Eq.  (68)  were  calculated  using 
either  Eq.  (1)  or  Eq.  (2).  Values  of  were  read  off  Figs.  (2),  (3) 
or  (4)  assuming  a  relative  velocity  between  the  fuel  drop  and  the  surrounding 
gas  of  twenty-five  percent  of  the  gas  velocity,  U^. 
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For  the  three  tubular  combustors  examined,  namely  the  J79-17A, 

J79-17C,  and  TF41 ,  values  of  Q'  of  0.99,  1.03  and  1.35  respectively 
provided  excellent  correlations  of  the  experimental  data,  as  demonstrated 
in  Figs.  70,  71  and  72.  These  values  of  Q’  correspond  to  liner 
pressure  loss  factors  for  these  combustors  of  14,  15  and  19  respectively, 
all  of  which  seem  quite  reasonable.  It  is  of  interest  to  note  that  the 
improvement  in  pattern  factor  with  increase  in  engine  power,  as  predicted 
by  Eq.  (68),  (due  to  reduction  in  evaporation  time),  is  fully  borne  out 
by  the  results  contained  in  Figs.  70,  71  and  72. 

The  influence  of  fuel  type  on  pattern  factor  is  manifested  through  the 
effects  of  mean  drop  size  (via  viscosity  and  surface  tension)  and  effective 
evaporation  constant  (via  T^)  on  droplet  evaporation  time.  Over  the 
range  of  fuels  examined  the  effect  of  fuel  type  on  pattern  factor  is 
relatively  small,  at  least  at  high  power  conditions  where  the  evaporation 
time  is  always  a  small  fraction  of  the  total  combustor  residence  time, 
regardless  of  fuel  type.  However,  if  measurements  of  pattern  factor  are 
conducted  at  operating  conditions  where  the  evaporation  time  constitutes  a 
significant  proportion  of  the  total  residence  time,  then  a  strong  effect 
of  fuel  type  on  pattern  factor  should  be  expected.  This,  in  fact,  was 
precisely  the  result  obtained  with  the  F101  combustor  when  the  effect  of 
fuel  type  on  pattern  factor  was  examined  at  various  simulated  engine 
operating  conditions,  but  using  air  supplied  at  normal  atmospheric 
pressure.  The  results  of  these  tests  are  shown  in  Fig.  73,  where  it  is 
of  interest  to  note  that  the  measured  values  of  pattern  factor  are  well 
correlated  by  Eq.  (68)  using  a  value  for  Q'  of  2.0.  Figure  73  demonstrates 
a  clear  effect  of  fuel  type  on  pattern  factor  but  it  would  be  erroneous 
to  assume  that  these  data  have  any  relevance  to  the  real  engine.  If 
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Comparison  of  Measured  and  Predicted  Values  of  Patte 
Factor  for  J79-17A  Combustor. 


Eq.  (68)  is  used  to  calculate  values  of  pattern  factor  at  actual  take-off 
conditions,  including  P^,  it  is  found  that  pattern  factor  is  virtually 
independent  of  fuel  type,  as  shown  in  Fig.  ’’4.  Thus,  from  a  practical 
viewpoint,  the  results  of  the  F101  tests  are  fully  consistent  with  those 
of  the  J79-17A,  J79-17C,  and  the  TF41 .  They  all  confirm  that,  at  the 
high  power  conditions  where  pattern  factor  is  most  important  to  engine 
durability,  variation  in  fuel  type  has  a  negligible  effect. 
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SECTION  XI 


DISCUSSION  AND  SUMMARY 


Analysis  of  the  key  processes  occurring  within  gas  turbine  combustors, 
along  with  examination  of  the  experimental  data  contained  in  references 
1  to  6,  shows  that  although  tne  impact  of  fuel  type  on  combustion  performance 
and  liner  durability  is  usually  small  in  comparison  with  the  effects  of 
liner  geometry  and  combustor  operating  conditions,  it  is  nevertheless 
of  sufficient  magnitude  to  warrant  serious  consideration.  For  some 
performance  parameters,  such  as  liner  wall  temperature  and  exhaust  smoke, 
it  is  found  that  fuel  chemistry  plays  an  important  role.  For  others, 
the  effects  of  fuel  type  are  manifested  through  the  physical  properties 
that  govern  atomization  quality  and  evaporation  rates. 

In  the  following  sections  the  effects  of  liner  size,  liner  pressure 
drop,  combustor  operating  conditions,  and  fuel  type  on  various  aspects  of 
combustion  performance  are  briefly  reviewed  in  turn. 


1 .  Combustion  Efficiency 

For  reaction-rate-controlled  systems 

P,1'75  Vc  exp(T3/300) 

* 

For  mixmg-controlled  systems 


(5) 


(8) 


For  evaporation-rate-controlled  systems 


*eff  ‘g  Vc 


(14) 
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From  analysis  of  the  experimental  data  contained  in  references  1  to  6  it 
was  found  that  combustion  efficiency  is  obtained  as  the  product  of  the 


0  efficiency,  n  ,  and  the  evaporation  efficiency,  nr  >  i.e. 

co  e 


nc  ”  nc  x  nc 
ce  ce 


where  nc  =  1  -  exp 


1.3 


-0.022  P3'-J  Vc  exp  (Tc/400) 


fc  mA 


(22) 


and  nc  =  1  -  exp 


■36  x  10J  P3  Vc  Xeff 


Tc  V  fcmA 


(23) 


Inspection  of  these  equations  reveals  that  combustion  efficiency  is 
enhanced  by  increases  in  P3,  T3  (via  Tc) ,  and  Vc>  and  by  reduction  in 
combustor  air  flow  rate.  There  is  a  very  slight  effect  of  fuel  chemistry 
on  nc.  via  Tc>  but  the  main  influence  of  fuel  type  on  combustion  efficiency 
stems  from  the  physical  properties  up>  op  and  T^n  that  govern  the  values 
of  DQ  and  Xeff  in  Eq.  (23). 


2.  Lean  Blowout 

Weak  extinction  values  of  fuel/air  ratio  are  obtained  as 


A'  f  m„ 
pz  A 


qLB0 


Vc  p31,3  exp(T3/300) 


V  LCVr 


-  x 


D  at  Tr 


D  at  277.5  K 


(33) 


In  this  equation  it  is  of  interest  to  note  that  the  dependence  of 
weak  extinction  limits  on  combustor  volume  and  operating  conditions  is 
identical  to  that  for  combustion  efficiency.  Also  in  common  with 
combustion  efficiency  is  the  slight  effect  of  fuel  chemistry  (via  LCVr) , 
whereas  physical  properties  are  important  due  to  their  influence  on  Dr 
and 

3.  Lean  Liqhtup 

The  equation  for  lean  lightup  fuel/air  ratio  is  virtually  identical 
to  that  for  lean  blowout  fuel/air  ratio,  except  for  a  slightly  stronger 
dependence  on  Pg.  We  have 

B'fp>A  Dr2  ~DatTF  "'2  .... 

qLL0  '  VCV5  exp  !T3/300>  X  V  LCVr  ”  1“  “  777'5KJ  <3  ’ 

Thus  the  observations  made  above  on  the  influence  of  combustor 
operating  conditions,  and  the  relative  importance  of  physical  and  chemical 
fuel  properties,  on  combustion  efficiency  and  lean  blowout  limits,  apply 
with  equal  force  to  ignition  performance. 

4.  Liner  Wall  Temperature 

The  most  important  factor  governing  liner  wall  temperature  is  the 
combustor  inlet  temperature,  T^.  Inlet  pressure  is  also  significant 
due  to  its  influence  on  the  concentration  of  soot  particles  in  the  flame, 
and  hence  on  the  magnitude  of  the  luminous  radiation  flux  to  the  liner 
wall.  At  max  power  conditions,  where  liner  wall  temperatures  are  of 
most  concern,  evaporation  rates  are  so  high  that  the  physical  properties 
of  the  fuel  appear  to  have  a  negligible  influence  on  T  .  Chemical 
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effects  are  also  quite  small,  as  shown  in  Figs.  44  to  47.  However,  even 
small  increases  in  maximum  values  of  liner  wall  temperature  can  seriously 
curtail  liner  life.  Thus,  for  the  range  of  fuels  covered  in  this 
investigation,  fuel  type  must  be  considered  of  significance  to  liner 
durability. 

In  the  calculation  of  liner  wall  temperatures  the  effect  of  fuel  type 
can  be  accommoda ted  quite  conveniently  by  introducing  the  fuel  hydrogen 
content  into  the  existing  equation  for  gas  emissivity.  This  approach 
leads  to  the  following  equation  for  eg. 

eg  =  1  -  exp  -97.44  P3(%H2)-2  (q  *b)0*5  Tg"1>5  (50) 

5.  NO  Emissions 

X _ 

It  is  found  that  N.0X  emissions  are  very  dependent  on  combustor 
operating  conditions,  and  also  on  the  size  of  the  combustion  zone  which 
governs  the  time  available  for  N0X  formation.  The  key  factor  controlling 
N0x  is  the  stoichiometric  flame  temperature  which,  in  turn,  is  almost 
solely  dependent  on  combustor  inlet  temperature.  As  far  as  fuel  type  is 
concerned,  physical  properties  are  of  little  consequence  except  at  low 
power  conditions  where  N0X  emissions  are  always  quite  small  due  to  the 
correspondingly  low  values  of  T  .  Fuel  chemistry  also  has  little 
influence  on  N0X  because  it  affects  only  slightly  the  values  of  bulk  gas 
temperature  (Tpz),  and  stoichiometric  flame  temperature  (Tst),  in  the 
following  equation  for  N0y. 

9  x  10'8  pj*25  V„  exp  (0.01  T  J 
N°x  = - - - - g/kg  (54) 

mA  Tpz 
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6.  CO  Emissions 

These  are  correlated  by  the  expression:  - 


36  "A  Tpz  exp  (-0-00345  TPZ) 
- H - • -  71 -  /T 

vc  -  o.55  5-—r~  &] 

[  Ppz  effj  P3 


It  is  again  observed  that  combustor  size  and  operating  conditions 
play  a  prominent  role  in  determining  the  level  of  CO  emissions.  Special 
importance  is  attached  to  inlet  temperature  and  primary-zone  fuel/air 
ratio,  due  to  their  combined  effect  in  resolving  the  primary-zone  temperature. 
As  in  the  case  of  NOx  emissions,  the  influence  of  fuel  chemistry  is  small 
and  is  manifested  through  slight  variations  in  Tpz  with  changes  in  lower 
calorific  value.  However,  since  CO  emissions  attain  their  maximum  values 
at  low  power  conditions,  where  a  sizeable  proportion  of  the  total 
residence  time  in  the  combustion  zone  is  occupied  by  evaporation  processes, 
the  influence  of  those  physical  properties  which  affect  evaporation  rates, 
namely  up,  op  and  Tbn,  becomes  important.  On  this  basis  it  would  be 
anticipated  that  fuels  of  high  viscosity  would  be  characterized  by 
slightly  higher  levels  of  CO  emissions,  and  the  experimental  data  generally 
confirm  this  expectation. 


7.  Smoke 


Of  all  the  parameters  studied,  smoke  emissions  is  the  one  that  is  most 
affected  by  changes  in  fuel  type.  The  physical  properties  of  the  fuel  are 
important  insofar  as  they  influence  the  mean  drop  size  in  the  spray  and 
the  penetration  of  the  spray  across  the  combustion  zone.  Spray  penetration 
is  of  considerable  importance  from  a  smoke  viewpoint  because  inadequate 
penetration  leads  to  enhanced  fuel  enrichment  of  the  soot-forming  regions 
just  downstream  of  the  fuel  injector.  Unfortunately,  present  knowledge  of 
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spray  characteristics  is  insufficient  to  allow  these  effects  to  be  described 
quantitatively.  However,  as  mean  drop  size  and  spray  penetration 
(especially  for  swirl  pressure  atomizers)  are  both  dependent  on  pressure, 
their  combined  effect  can  be  embodied  into  a  pressure  term  such  as  that 
contained  in  the  following  equation  for  soot  concentration  in  the 
combustor  exhaust  gas. 


X  *  v . 

fpz  Tpz 


A  -  exp  (0.0014  TS2) 


Aromatic  Content  +  B 


(60) 


The  above  equation  shows  that  combustion  pressure  and  primary-zone 
fuel/air  ratio  are  the  main  factors  governing  the  output  of  smoke. 

However,  fuel  type  also  has  a  strong  influence,  as  evidenced  by  the 
relationships  between  soot  concentration  and  aromatics  content,  shown 
plotted  in  Figs.  64  to  67.  Although  total  aromatics  content  was  selected 
as  the  best  indicator  of  a  fuels'  propensity  for  soot  formation,  the  form 
in  which  it  is  included  in  Eq.  (60)  is  unsatisfactory,  since  different 
combustors  require  different  values  of  B.  More  work  is  needed  to  clarify 
the  role  of  fuel  composition  and  combustion  chemistry  on  soot  formation 
and  smoke.  Carefully  controlled  experiments,  conducted  at  high  pressures 
on  selected  fuels,  incorporating  accurate  measurements  of  spray  character¬ 
istics,  should  provide  a  more  definite  indication  of  the  relative  importance 
of  hydrogen  content  and  aromatics  type  than  can  be  deduced  from  the  existing 
data. 


8.  Pattern  Factor 

This  is  described  with  good  accuracy  by  the  following  equation 


T 


max 


T 


4 


-  exp  -  Q' 


r 

|ll  • 

0.33  mfl  o/ 

^ref 

pg  AL°LXeff 

(68) 
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where  appropriate  values  of  Q1  are  0.70  and  0.50  for  tubo-annular  and 
annular  combustors  respectively.  The  above  equation  shows  that  the  two 
main  parameters  controlling  pattern  factor  are  the  pressure  drop  across 
the  liner  and  the  liner  L/D  ratio.  It  also  accounts  for  the  influence  of 
evaporation  time  in  reducing  the  time  available  for  mixing  within  the 
liner.  At  the  high  pressure  conditions  where  pattern  factor  is  of  most 
concern,  the  evaporation  time  is  always  quite  short  in  comparison  to  the 
total  residence  time  of  the  combustor,  and  so  the  dependence  of  pattern 
factor  on  fuel  type  is  fairly  small,  as  illustrated  in  Fig.  74. 

With  reduction  in  engine  power  the  evaporation  time  increases  due 
to  increase  in  0Q  and  reduction  in  This  produces  a  deterioration 

in  pattern  factor  as  indicated  by  Eq.  (68)  and  also  by  Figs.  70  to  73, 
which  demonstrate  for  all  engines  that  pattern  factor  at  idle  is  distinctly 
worse  than  at  take-off.  These  considerations  highlight  the  importance  of 
measuring  pattern  factor  only  at  the  correct  combustor  inlet  conditions 
of  mA,  P3>  T3  and  qQV,  corresponding  to  engine  operation  at  max  power. 

Tests  carried  out  at  simulated  conditions  at  lower  pressure  levels  give 
misleading  results,  as  shown  in  Fig.  73.  First,  they  yield  values  that 
are  overly  pessimistic  and,  second,  they  show  a  dependence  of  pattern  factor 

on  fuel  type  which  greatly  exaggerates  the  dependence  actually  observed 
at  high  pressure. 


SECTION  XII 
CONCLUSIONS 

1.  The  most  important  factors  governing  liner  durability  and  combustion 
performance  are  liner  size,  liner  pressure  drop  and  combustor  operating 
conditions,  with  fuel  effects  playing  a  secondary  role.  However, 
since  in  modem  high  pressure  ratio  engines  the  combustor  is  called 
upon  to  perform  satisfactorily  for  long  periods  at  extreme  conditions 
on  current  fuels,  it  follows  that  any  factor,  however  secondary,  that 
creates  a  more  adverse  combustor  environment,  will  have  a  large  and 
disproportionate  effect  on  combustion  performance  and  liner  durability. 

2.  Analysis  of  the  experimental  data,  which  cover  a  range  of  fuel  types 
from  JP4  to  0F2,  shows  that  fuel  chemistry,  as  Indicated  by  hydrogen 
content  and/or  aromatics  content,  has  a  significant  effect  on  flame 
radiation,  liner  wall  temperature  and  smoke  emissions. 

3.  The  influence  of  fuel  chemistry  on  Ignition  performance,  weak  extinction 
limits,  combustion  efficiency,  and  CO  and  N0X  emissions,  is  quite  small, 
and  stems  from  the  effects  of  slight  variations  in  lower  calorific  value 
on  combustion  temperature. 

4.  The  physical  properties  that  govern  atomization  quality  and  evaporation 
rates  affect  lightup  characteristics,  weak  extinction  limits,  combustion 
efficiency,  and  CO  emissions.  Other  important  performance  parameters, 
such  as  N0X  emissions,  smoke  emissions  and  liner  wall  temperature,  are 
sensibly  Independent  of  physical  properties  over  the  range  of  fuels 
studied. 


Preceding  Page  Blank 
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5.  Fuel  chemistry  has  no  direct  influence  on  pattern  factor.  However, 
physical  properties  have  an  effect  that  is  appreciable  at  low  power 
conditions  but  which  diminishes  in  importance  with  increase  in  engine 
power,  becoming  very  small  at  the  highest  power  setting,  where  the 
effect  of  pattern  factor  on  vane  life  is  most  significant. 


SECTION  XIII 


RECOMMENDATIONS 

1.  Additional  work  is  recommended  to  extend  the  scope  of  the  present 

report  to  include  correlations  for  the  prediction  of  unburned  hydrocarbons 
(LHC).  During  this  program  the  data  on  unburned  hydrocarbons  contained 
in  references  1  to  6  were  examined  to  the  extent  that  it  was  found 
feasible  to  obtain  correlations,  although  not  to  the  same  degree  of 
accuracy  as  for  CO  and  N0y.  However,  lack  of  time  precluded  the 
establishment  of  a  satisfactory  correlation  for  UHC. 

2.  The  correlations  presented  in  this  report  should  be  re-examined  with 
a  view  to  further  improvement  and  refinement.  For  example,  the 
present  soot  correlation  is  based  on  the  total  aromatic  content  of 
the  fuel.  However,  it  is  quite  possible  that  not  all  of  the  various 
types  of  aromatics  (1-ring,  2-ring,  etc.)  are  equally  effective  in 
producing  soot;  and  thus  a  better  correlation  of  smoke  emissions 
might  be  achieved  by  neglecting  certain  aromatic  constituents  of 
the  fuel.  Also  worthy  of  further  investigation  is  the  use  of 
hydrogen  content,  instead  of  percentage  aromatics,  as  an  indicator 
of  propensity  to  soot  formation. 

3.  It  is  strongly  recommended  that  efforts  be  made  to  measure  the  mean 
drop  sizes  (SMD's)  produced  by  the  fuel  nozzles  actually  used  in 

the  fuels  research  programs  described  in  references  1  to  6.  Accurate 
knowledge  of  drop  sizes  over  the  entire  engine  operating  range  would 
not  only  allow  the  fairly  complete  correlations  achieved  for  the 
J79-17A,  J79-17C,  F101,  and  TF41  to  be  applied  to  the  TF39,  J85,  TF33 
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and  FI 00,  but  would  also  enable  these  correlations  to  be  rewritten 
in  a  more  basic  and  general  form  that  would  facilitate  their  application 
tc  all  other  gas  turbine  combustors. 
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